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Abstract 

An extensive literature exists describing various techniques for the evaluation of Hankel deter- 
minants. The prevailing methods such as Dodgson condensation, continued fraction expansion, LU 
decomposition, all produce product formulas when they are applicable. We mention the classic case 
of the Hankel determinants with binomial entries ( 3fc ^ 2 ) and those with entries f 3 ,?) ; both of these 
classes of Hankel determinants have product form evaluations. The intermediate case, ( 3fc fc l_1 ) has 
not been evaluated. There is a good reason for this: these latter determinants do not have product 
form evaluations. 

In this paper we evaluate the Hankel determinant of ( 3k ^ 1 ) ■ The evaluation is a sum of a small 
number of products, an almost product. The method actually provides more, and as applications, 
we present the salient points for the evaluation of a number of other Hankel determinants with 
polynomial entries, along with product and almost product form evaluations at special points. 



1 Introduction 

A determinant 

H n = det[dij]o<i,j<n 

whose entries satisfy 

for some sequence {ak}k>o is said to be a Hankel determinant. Thus H n is the determinant of a special 
type of (re + 1) x (n + 1) symmetric matrix. 
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In various cases of Hankel determinant evaluations, special techniques such as Dodgson conden- 
sation, continued fraction expansion, and LU decomposition are applicable. These methods provide 
product formulas for a large class of Hankel determinants. A modern treatment of the theory of deter- 
minant evaluation including Hankel determinants as well as a substantial bibliography can be found 
in Krattenthaler [5j [6] . 

The product form determinants of special note are those whose factors have some particular at- 
traction. Factorials and other familiar combinatorial entities that appear as factors have an especially 
pleasing quality, and we find an extensive literature devoted to the evaluation of classes of Hankel 
determinants as such products. 

Several classical Hankel determinants involve entries that are binomial coefficients or expressions 
closely related to binomial coefficients. Perhaps the most well-known of these is where a& = ( 2fc fc l_1 ), 
and a/; = 2k+i C k k~ l ) ' f° r wn i cn Hn = 1 for all n. 

It is also known that binomial entries such as 



yield product evaluations for the corresponding H n . In fact, product formulas have been shown to 
exist for a host of other cases (see Gessel and Xin [2]), and we only mention 



as representatives. 

However, within the restricted class of Hankel determinants defined by the binomial coefficients 



parametrized by a pair of integers (3 > and a, it is a rare phenomenon that the determinant 
evaluations are in product form. An extensive check of Hankel determinants of sequences in the 
form ([2]) suggests that there is no product formula for H n in general for such binomial sequences. In 
fact, it would seem that the instances for which H n has a product form can be enumerated in full: 







(2) 



(i) (3 = 1, a arbitrary, 



(ii) (3 = 2, a = 0,1,2,3,4, 



(iii) = 3, a = 0,2. 
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All the other cases are likely not products, but in any event, the question remains open: are 
evaluations possible in these cases? 

Let H^' a ^ denote the (n + 1) x (n + 1) Hankel determinant with entries ajjf' 1 ^ as defined in ([2]). 
We will also use the term ((3, a) -case to refer to the evaluation of H^ ,a \ 

Numerical data indicates the intriguing possibility that the H^'°^ might be evaluated as a sum of 
a small number of products, where "small" would mean 0(n d ) summands for some fixed d = d(f3, a). 
We refer to such an evaluation as an almost product. 

The evidence of an almost product evaluation of H^ 3 ' ^ is most pronounced for (3 = 3, and we 
begin with some sample data. 

For the (3, 2)-case the Hankel determinants evaluate to 



H (3,2) = 2 2 ■ 3 ■ 7 3 ■ 37 • 41 2 ■ 43 3 ■ 47 3 ■ 53 2 -59-61 

o-(3,2) 
"20 

H fo 2) = 2 10 ■ 5 12 ■ ll 9 ■ 13 3 • 41 3 ■ 43 3 ■ 97 ■ 101 2 ■ 103 3 • 107 4 ■ 109 s ■ 113 6 • 127*0 ■ 131 9 ■ 137 s ■ 139 s ■ 149 6 • 
151 6 • 157 s • 163 4 • 167 3 • 173 2 ■ 179 • 181 

The small prime factors are indicative of the fact that there is an underlying product formula. In fact 
for the (3, 2)-case, the Hankel determinant is explicitly given by 

H (3,2) = n ^+m^ + iy. m 

^ 2(4i + 2)!(4i + 3)! 1 ; 

We mention also the (3, 0)-case for which the Hankel determinant also has small prime factors, 
and possesses the product evaluation 

m " 3(3i + l)(6i)!(2i)! 

l[ (4»)l(4» + l)! U 

For the (3, l)-case, we get the following intriguing evaluations: 

H[l' 1} = 2 2 • 7 2 • 37 • 41 2 ■ 43 3 ■ 47 2 • 53 • 41740796329 

i^ 3 ' 1 ' = 3 s ■ 29 • 67 • 71 2 • 73 3 • 79 5 • 83 6 ■ 89 5 ■ 97 4 • 101 3 ■ 103 3 ■ 107 2 • 109 2 • 113 • 631 ■ 548377971864917477341 
= 2 10 ■ 5 10 ■ ll 9 ■ 13 2 ■ 41 3 • 43 2 ■ 97 • 101 2 ■ 103 3 ■ 107 4 • 109 5 • 113 6 ■ 127 9 • 131 s • 137 7 ■ 139 7 ■ 149 5 • 
151 5 ■ 157 4 • 163 3 • 167 2 • 173 • 569 • 920397320923 • 56029201596264233691799 

(3 1) 

The existence of large primes in the factorizations indicates that Hn ' does not have a product form 
evaluation. However if we write H n = P n Q n where P n is the product of the small primes and Qn 
is the product of the large primes left over, then it appears that the estimates logP n = Q(n 2 ) and 
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logQn = 0(n) hold. This suggests that these Hankel determinants can be represented as a sum of 
0(n) number of products, all of which have very similar representations. 

(3 1) 

The purpose of this paper is to provide a method that evaluates Hn ' and a number of other 

(3 1) 

Hankel determinants as almost products. For H n = H n , we obtain 

" (6i-3)!(3i + 2)!(2i-l)! " n!(3n + i + 2)!(-6) i 
n i 1 = i 1 (4i-l)!(4i + l)!(3i-2)!^(3n + 2)!(n-z)!(2i + l)!' {> 

or alternately 

" (6i + 4)!(2i + 1)! " n!(4n + 3)!!(3n + i + 2)! 
" ~ AX 2 (4i + 2)!(4i + 3)! ^ (3n + 2)!*!(n - i)\(in + 2i + 3)!! ' ^ ' 

each as a sum of n + 1 products. The method actually evaluates more, and we describe now the general 

situation in the (3, l)-case, which consists of three basic ingredients: 

(I) Replace a k with polynomials 

»,(x)=4 3 ' 1, w=i:( 3 \ + _ 1 ;i-™ m 

so that cifc(x) is a monic polynomial of degree k with = afc(O). 
(II) Show that the dk{x) satisfy certain differential-convolution equations. 

(Ill) Show that the resulting determinants H n (x) themselves satisfy certain differential equations. 



(3 1) 

The (n + 1) x (n + 1) Hankel determinant H n (x) = H n ' (x) is then expressed as the power series 
solution of the differential equation in (III), and we give it here as it is stated as Theorem [2j 

H( \ -( IV* TT (6i-3)!(3i + 2)!(2i-l)! " n!(3n + z + 2)!2^ (x - 3)^ 
n (4i - 1)!(4» + l)!(3i - 2)! ^ (3ra + 2)!(n - i)!(2i + 1)!' { > 

An alternate expression for this evaluation appears in Theorem 1101 in Section [HI 

Similarly, steps (I), (II), (III), mutatis mutandis, yield (Theorem [5]) : 

H^(x) - (-IT (2n + 1)T (n + ^ ( "" 2)t 
H n {x)-{ L) Vn + L)^ {n _ i)]{ 2i + 1)1 ' 

where 

k 



a k {x) = a^\x)=Y, [ U m - 0) 

These polynomial families have a number of interesting properties that we briefly discuss. For 

(3 1) 

example, the polynomials H n ' (x) satisfy a three-term recursion, their roots are real, and interlace. 
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Furthermore the specializations at x = 3, §, | all have product evaluations. The polynomials Hn ' (x) 
form an orthogonal family. A few other classes of Hankel determinants can be evaluated in almost 
product form by simple transformations of these polynomials (e.g., Example [HJ in Section [HJ and 
Corollaries [7J and [HJ in Section 16. 2\) . 

Returning briefly to the general case of the determinants H^' a \ we remark that there are a few 
more cases which fall under the method described above. These would also include the same three 
ingredients: 

(I) Replace with polynomials 



a k {x) = 




so that a^(x) is a monic polynomial of degree k with a k = Ofc(0). 
(II) Show that the a,k(x) satisfy certain differential-convolution equations. 

(Ill) Show that the resulting determinants H n (x) themselves satisfy certain differential equations. 

The (3, l)-case and the (2, l)-case are governed by second order differential equations, but even 
these cases already present considerable technical problems to overcome. We mention some further 
difficulties that arise in the consideration of other (/3, a)-cases in Section [7J A number of additional 
almost product evaluations of Hankel determinants are given in Section [HJ For these additional results 
given as Theorems IS] \7\ [HI l9l [TOl and special product evaluations that appear in (|167p . (| 168|) and (|17ip . 
we provide the necessary identities for proving the differential equations, mimicking the proofs we 
present for Hn'^ix) and Hn'^^x). 

Finally, a remarkable property of these (n + 1) x (n + 1) Hankel determinants H^' a \x) is that the 
degree of the polynomial H n (x) is only n, indicating an extraordinary amount of cancellation in the 
expansion of the determinant. The unusual degree of cancellation is a basic property of a large class 
of Hankel determinants with polynomial entries. This class contains the Hankel determinants defined 
by polynomials in (fTUj) that we consider. The degree result is of independent interest, and we include 
an exact statement and a proof of it in Appendix III. 

We would like to remark that the differential-convolution equations (II) used in this paper are 
reminiscent of the equations that arise in the study of the Painleve II equation and the Toda lattice 

mm- 
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2 The (3, l)-case 

2.1 Differential-convolution equations 

In the proof of the (3, l)-case, we denote the polynomials a k ' (x) by a k , H n ' (x) by H n , and the 
differentiation operator by d x . 

We need two identities given below in Lemmas [1] and [2j The first is a differential-convolution 
equation. The second identity involves convolutions and a k but no derivatives. The proofs are given 
in Appendix II. 

Lemma 1 Let the polynomials a k = a k ' (x) be as defined in §7ty. Then 

(x - 3) (2a; - 3)(4x - 3)d x a k = 2{2k + 3)a fc+ i - (8x 2 - 18x + 27k + 36)a fe (11) 

+ 4(2x 2 - 6x + 3)c fc - 27(2x 2 - 6x + 3)c fc _i 

where 

k 

Ck = Cfc(x) = ^ a m (x)a fc _ m (x), (c_i = 0). (12) 

m=0 

Lemma 2 VFii/i = ' (x) and Ck = Ck(x) as in LemmaUl we have 

4(2k + 5)(x - l)a fc+2 - (2(16x 3 - 72x 2 + 135x - 81)fc + 2(24x 3 - 92x 2 + 180x - 117))a fc+ i 
+ (27(2x - 3) 3 A: + 54(2x - 3)(2x 2 - 4x + 3))o fc + 8(x - l)(2x 2 - 6x + 3)c k+1 
+ 2(8x 4 - 114x 3 + 324x 2 - 297x + 81)c fc - 27x(2x - 3)(2x 2 - 12x + 9)c fc _i = 

Lemmas [1] and [2] will be needed for the proof of the differential equation satisfied by the determi- 
nants H n . This differential equation is given below in Theorem [TJ 

In addition to the first two identities in Lemma [1] and [21 a much more complicated third identity 
involving the a k is also needed for the proof of this differential equation. This will emerge in the course 
of the proof of (|T4"|) . 

Theorem 1 Let the polynomials a k = a k ' (x) be as in §7ty and define the (n + 1) x (n + 1) Hankel 
matrix by 

A n = A n (x) = [ai+j(x)]o<i,j<n 

Then 

H n = H n 3 ' 1 \x)=detA n (x) (13) 

satisfies the differential equation 

(x - l)(x - 3)d 2 x y + (2(n + 2)(x - 3) + 3) d x y - 3n(n + l)y = 0. (14) 
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Proof This is easy to check for n = 0, 1. Henceforth we assume that n > 2. 

We will find expression for the first derivative d x H n Section 12.21 and the second derivative d x H n 
Section [2.31 This is followed by Section [3] on specializations of x: product formulas for H n (3), H n (^) 
and H n (j) are given as three corollaries in Sections 13.11 13.21 and 13.31 The specializations make use 
of a Dodgson-like expansion result we prove as Proposition [1] at the start of Section O The proof of 
Theorem [T] continues in Section [3] where we put together the expressions obtained for the derivatives 
and the third identity mentioned to prove (|14p . 

2.2 Calculating the First Derivative 

The first step is to find a simple form for the derivative of H n . We begin with the expression 

d x H n = Ti{A- l d x A n )H n (15) 
for the derivative of a determinant, where 

d x A n = d x A n (x) = [d x aj+j (x)]o<i,j<n- 
Referring to Lemma [U we write 

(x - 3)(2x - 3)(4x - 3)Tr(A- 1 4.A„) 

as 



2Tr(^ 1 [(2(i+j) + 3)a i+j+l]o<i,j<n) (16) 

+Tr(^ 1 [-(8x 2 - I8x + 27(i + j) + 36)a i+i ]o<i,i<n) (17) 

+4(2x 2 - 6x + 3)Tr(A- x [ci +j ] < id < n ) (18) 

-27(2x 2 - 6x + 3)Tr(A" 1 [c i+J _i] <i j<„) (19) 

where the convolutions are defined in (|12p . 



We render each of these four expressions f)16[> -(|19 |) in a simple form, and then combine them all 
into an expression for the derivative in (|15p . After that is done, we go through a similar computation 
for the second derivative of H n , where we use the recursion in Lemma [2] for the simplifications. The 
differential equation will follow from a third identity, the proof of which makes up the bulk of the work 
for the rest of the argument. 
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We begin with the trace term (| 16[) : Let I denote the identity matrix of relevant dimension and 
define the two matrices 



B n = B n (x) 



L, 



Then 
2Tr 



(A' 1 [(2(i+j)+3)a l+j+1 (x)} 

0<i,j<n 

Now define ctq, oi, . . . , a n and K n as follows 



a\ a,2 ... a n+ i 

02 «3 ... a n+ 2 

O-n+1 &n+2 ■ ■ ■ «2n+l 


1 

2 

n 



2Tr ( A- 1 ((2L n + h)B n + 5„(2L n + 



and 



K n = det 



0"o 










a n +2 


_ Cn _ 






a 


ai 


a n _i 


a i 


a 2 





(20) 



O n -l «n • • • 02n-2 «2n 
On On+1 • • • 02n-l a2n+l 



(21) 



By Cramer's rule we have 



0"„ 



(22) 
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Therefore 



A n B n 



and 



A- 1 B n (2L n + -I) 



... ao 

1 o-i 

1 (J 2 

cr n _i 

1 CJ n 



... (2n + 3/2)a 

3/2 (2n + 3/2)cri 

7/2 (2n + 3/2)cr 2 

(2n + 3/2)a n _i 

(4n-l)/2 (2n + 3/2)o-„ 



Since 



we can write 



(2L n + -I)B n A~ l = {A~ l B n (2L n + -I)) T 



2Tr ( A^((2L n + h)B n + 5„,(2L n + ^7)) 

2Tr (j^\2L n + + 2Tr ^- 1 S„,(2L n + |j) 

2Tr {{2L n + -^T^" 1 ) + 2Tr (A- l B n (2L n + ^7) 



4Tr ^- 1 J B n (2L n + -I) 
3 

4(2n + -)a n 



2(4n + 3)a-„ 



The last expression gives 



2Ti(A- l [{2{i+j) + 3)a t+j+1 ] < Lj < n ) = 2(4n + 3)^. 



for the desired form of the first term (|16p . It is useful to record in passing that 



The identity (I25h will be useful later when we calculate the derivative of K n . 



Now we consider the second trace term (|17p . The calculation of this term is done in the same 
manner as the first term but the evaluation is somewhat simpler. We get the expression: 



Tr^A^l-fix 2 - 18x + 27(i + j) + 35)a 



i+j\0<i,j<n) 



-(n + l)(8x 2 - 18x + 27n + 36). (26) 



The final two terms (|18p and (|19p require a new technique. We will use ideas from [H [3] where an 
identity similar to the following is used: 



-i+j\0<i,j< 



n — E n A n + A n E^ 



(27) 



where 



E n = EJx) 



a /2 
ai/2 a 
CL2/2 a\ a 



(28) 



On/2 a n -i a n -2 ao 
Note that the first column of E n is divided by two. This allows for the immediate computation 

Tr^n 1 [ci+i] < w <„) = ( 2n + ijoo- 

So the third term (|18p is 

4(2x 2 - 6x + 3)Tr(A^ 1 [c i+J ] <i, j <„) = 4(2n + l)(2x 2 - 6x + 3) 



(29) 



The term (11911 that involves 



l c i+j-l\o<ij< n 



case is similar. We have the equation 



[ c i+j-l\o<i,j<n ~ FnA n + A n F n 



(30) 



where 



En — E n (a. - J 





a 

ai ao 

On-l «n~2 • • • «0 



(31) 



The identity (|30|) leads to the computation 

Tr(^' 1 [ci +3 -_ 1 ]o<ij<n) = 0, 
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qo that term (fl~9|) evaluates to zero: 

- 27(2x 2 - 6x + 3)TV(^- 1 [c, +j _i] < lJ <n) 







(32) 



Adding the expressions (f24"|) , ([26]) , ([29]) , ([32]) and multiplying through by i/ n we obtain the following 
expression for the first derivative 

Lemma 3 

(x - 3)(2x - 3) {Ax - 3)d x H n = (8nx 2 - 6(5n + l)x - 3(9n 2 + 13n + 8)) H n + 2(4n + 3)K n (33) 

We next state and prove a lemma which is preparatory to the calculation of the second derivative 
of H n . First define two new determinants as follows: 



M n = M n (x) = det 



N n = N n (x) = det 



a 


ai 


«n-2 


On+1 


a„ 




a 2 


«n-l 


O n +2 


O-n+l 




a n 


• • A2n-3 


«2n 


a 2n -i 


a n 




• • Gt2n-2 


«2n+l 


0-2n 


«o 


a\ 


O-n-2 




O-n+2 


«i 


a 2 






a n +3 




a n 


• «2n-3 


^2n-2 


02n+l 


a n 


a n +i ■ 


• «2n-2 


d2n-l 


0.2n+2 



(34) 



(35) 



where = a^(x) are the polynomials defined in (|7j). Then there is a linear relationship between the 
four determinants H n , K n , M n , and N n as stated in the following lemma. 



Lemma 4 



4(4n + 5) (a - l)JV n + 4(4n + l)(x - l)M n 
+ (-16(4n + l)x 3 + 8(36n + 7)x 2 - 108(5n + 2)x + 6(54n + 31)) K n 
+ ((64n + 16)x 4 + (216n 2 - 24n - 12)x 3 - (972n 2 + 800n + 108)x 2 

+(1458n 2 + 1770n + 378)x - 729n 2 - 1083n - 324) H n = 



(36) 



Proof of the Lemma: First we make use of the recursion in in Lemma [2] for each index k. Putting 
these all together in matrix form, we apply the operator 

Tr^- 1 *) 
11 



to obtain the trace identity 

Tr(A i ; 1 [4(2(z + j) + 5)(x - l)a i+j+2 ] < itj < n ) (37) 
+ Tr(4; 1 [(-2(16x 3 - 72a; 2 + 135.x - 81) (i + j) - 2(2Ax 3 - 92x 2 + 180x - 117))a l+i+ i] < lJ <„) (38) 

+ Tr(A- J [(27(2a; - 3) 3 ^ + j) + 54(2a; - 3){2x 2 - Ax + 3))a i+j ]o< iJ <„) (39) 

+ Tr(A- 1 [8(a; - l)(2x 2 - 6x + 3)c 4+ , +1 ] < lJ <„) (40) 

+ Tr^ 1 [2(8a; 4 - 114a; 3 + 324a; 2 - 297x + 81)c. i+j ] < 4 d < n ) (41) 

+ Tr(A,; 1 [-27x(2a; - 3) (2a; 2 - 12a + 9)c l+J _i] < lJ <„) - (42) 



Each of these six traces (|37p -(|42 p is calculated in a similar manner as was done above in the 
calculation of the first derivative of H n . The first and fourth trace will involve a small extension to 
what was used above. We will start with the computation of (|37p . As before we write this as 



4(s - l)Tr(^ 1 ((2L n + -I)[a i+j+2 }o<i,j<n + [a. l+j+2 }o<i,j<n{2L n + -/))). (43) 



Using the fact that 



4{x - l)Tt{A~ 1 {2L n + ll)[a i+j+2 ]o<i,j<n) 

= 4(x - l)Tr((2L n + |l)[a i+J - +2 ]o<i ! j<n^n 1 ) 



and 



(2L n + -I) [a i+ j + 2]o<i,j<nA n 1 j — A n 1 [ai + j +2 ]o<i,j<n(^L n + -I) 



we see that we can write (1431) as 



x - l)Ti{A n [a i+ j + 2]o<i,j<n(2L n + -J)). 



(44) 



We will obtain a representation of the trace above in terms of H n ,M n , and N n . Introduce to, n, . .., 
T n as follows: 



TO 






Tl 


— A n 


O ra +3 
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As before we observe that 







1 
1 



0"3 



70 

n 

73 



cr n _2 r n _ 2 
'"• (T n -1 T n _x 

1 cr« TV,, 



and therefore 



expands to 






5/2 



Thus (gU) simplifies to 
By Cramer's rule we have 



y ^n 1 [ a «+j+2]o<i,j<n(2£ n + -I) 



(An + l)a /2 (An + 5)r /2 

(4n + l)o-i/2 (4n + 5)Ti/2 

(4n + l)(r 2 /2 (4n + 5)r 2 /2 

(4n + IK-2/2 (4n + 5)r n _ 2 /2 

(4n + lK_i/2 (4n + 5)r n _!/2 

(4n - 3)/2 (4n + l)cr n /2 (4n + 5)r n /2 

4(x - l)((4n + l)o- n _i + (4n + 5)r„) 



c n -i 



so, in summary, the first trace (|37p evaluates to 



4(s - l)((4n + 1)^ + (An + 5)^). 



(45) 
(46) 

(47) 



Now we consider the second trace (j38|) . The evaluation of this trace proceeds exactly as the evaluation 
of the trace (fT6|) in the computation of the derivative of H n . ([38]) evaluates to 



2 (2(16x 3 - 72x 2 + 135x - 81)n + (24x 3 - 92x 2 + 180x - 117) 
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Kn 



(48) 



Similarly, the third trace (|39|) evaluates to 

(n + 1) (27 \2x - 3) 3 n + 54(2x - 3)(2x 2 - 4x + 3)) . (49) 

The next trace (|40p requires that we define a matrix G n in a similar manner to E n and F n . 
Specifically we define 

ai/2 cio/2 
02/2 ai/2 a 
03/2 02/2 a\ a 

a n /2 a n _i/2 a n _2 • • • 
On+i/2 a ra /2 a n _i . . . 



G n — G n (x) 



«o 

01 



With this definition of G n we can write 



[ c i+j+l]o<i,j<n — G n A n + A n G^ + do 



... 

O n +l O-n+2 ■ ■ ■ d2n+l 



+ a 



... a n+ i 
: a n+2 

... a 2n +i 



We have 



.4 



-1 



... a n+1 
: a n+2 



... a 2n +i 
where a,- is as defined in (l20l). It follows that 



... cr 

: o\ 

... a n 



( 



Tr 



A., 



... a n+ i 
■ a n+2 



... a 2r 



+1 J/ 



0"r).. 



The term 



Tr 



a: 



... 

O-n+l On+2 • • • d2n+l 
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comes out to a„ also because 



a: 



-i 



-i 



J 



a n+ i 
a n+2 







a 2 



n+l 



... 

On+1 a>n+2 ■ ■ ■ «2n+l 

Putting this together we see that the fourth trace (|40p evaluates to 

f ) ■ (50> 

The evaluation of the fifth and sixth trace (|4ip and ()42j) are done in exactly the same manner as 
the traces (1181) and (1191). For the fifth trace (1411) we obtain 



16(x - l)(2x iS - 6x + 3) n{x + 4) + 



(51) 



Tr(^" 1 [2(8x 4 - 114x 3 +324x 2 - 297x + %l)ci+j]Q<i,j<n) 

= 2(2n + l)(8x 4 - 114x 3 + 324x 2 - 297x + 81) 

and the sixth trace (|42p evaluates to zero just as the trace in Q19|) : 

Tr(^~ 1 [-27x(2x - 3)(2x 2 - 12x + 9)c i+J _i] <i, j <n) = . (52) 
Adding the expression we have found for the six traces in (|47 p -(|52 p we obtain the identity in (|36|) . 

2.3 Calculating the Second Derivative 

We are now ready to calculate the second derivative of H n . We begin with equation (I33p for the 
derivative. The first step is to replace K n in equation (133p with the representation of K n as a trace 
from (|25p . Inserting this expression into equation (|33p we get 



(x-3)(2x-3)(4x-3)4^n = (8nx 2 -6(5n + l)x-3(9n 2 + 13n + 8))H n + 2(4n + 3)Tv(A- 1 B n )H n (53) 

In order to obtain the second derivative of H n , we differentiate ([53]) . We obtain 

(x - 3)(2x - 3)(4x - 3)d 2 x H n + 
(8(n + 3)x 2 - 6(5n - 13)x - 3(9n 2 + 13n + 29))d x H n + 

2(8nx - 15n - 3)H n 
= 2(4n + 3)(d x Tr(A- l B n ))H n + 2(4n + 3)Tr(A- 1 B n )d x H n . 



Multiply both sides of the equation by 



(x-3)(2x-3)(4x-3) 
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The second trace term on the right hand side now becomes 

2(4n + 3)(x - 3)(2z - 3) (4s - 3)Tr (A' 1 B n )d x H n (54) 
Using the identities (|33|) and (|25p . (|54|) can be written in terms of H n and -ftT n as follows: as 

2(4n + 3) ^(8ns 2 - 6(5n + l)x - 3(9n 2 + 13n + 8))K n + 2(4n + 3)^-^J . 
Substituting back, we get 

(x - 3) 2 (2x - 3) 2 (4s - 3) 2 d 2 x H n - 
(x - 3) (2s - 3) (4s - 3)(8(n - 3)s 2 - 6(5n - 13)s - 3(9n 2 + 13n + 29))d x H n - 

2(x - 3) (2s - 3) (4s - 3)(8ns - 15n - 3)F n - 
2(4n + 3)(8ns 2 - 6(5n + l)s - 3(9n 2 + 13n + 8))K n - (55) 

4(4n + 3) 2 f- 



2(4n + 3)(s-3)(2s-3)(4s-3)(4Tr(yl~ 1 J B n )JF n = 
We will now focus on simplification of the derivative of the trace 

(s - 3) (2s - 3) (4s - 3)(d x Tr(^ 1 S n )) 
which is a factor of the last term on the left of equation (|55j) . We use the fact that 

dx-A n = —A n (d x A n )A n 

and write 

d^iA^Bn) = TiiA^dxBn) - i:T{A-\d x A n )A~ l B n ) (56) 

Using equation (fTTj) from Lemma [fl and multiplying equation ([5"o]) by (s — 3) (2s — 3) (4s — 3), we 
can write 

(s - 3) (2s - 3) (4s - 3)TT{A- x d x B n ) 



MS 



Tr(A- 1 [2(2(i + j + l)+3)a i+ j+2]o<i,j<n) (57) 

+ Tr(A- l [-(8x 2 - 18s + 27{i + j + 1) + 36)a i+J+1 ] < itj < n ) (58) 

+ Tr(A- x [4(2s 2 - 6s + 3)c J+i+1 ] < i j<„) (59) 

+ Tr^" 1 [-27(2s 2 - 6s + 3)c i+j } <i tj < n ) (60) 

16 



and we can write 

(x - 3)(2x - 3) (4s - 3)Tr(A-\d x A n )A- l B r , 



as 



Tr(^7 1 [2(2(f +j) + 3)a i+J+1 ] < ij <„^- 1 J B„) (61) 

+ Tv(A^[-(8x 2 - 18z + 27(* + i) + 36)a i+J ]o<ij< n ^n lB n) (62) 

+ Tr(^" 1 [4(2x 2 - 6x + 3)c J+i ]o< 4 tj <nA~ l B n ) (63) 

+ Tr(^~ 1 [-27(2x 2 - 6x + 3)c l+J _ 1 ] < l jXn^^J (64) 



We have already evaluated traces that are similar to the first four terms (|57p - (|60p in the calculation 
of the trace terms in (U6l)-(fT9l), In this case we obtain 



Tr(A- 1 [2(2( l + . 7 + l) + 3K +J+2 ]o<., J <„) = 2((4n + 1)^ + {An + 5& (65) 

Tr(A j ; 1 [-(8x 2 -18x + 27(i + j + l)+36)a, +J+ i]o< l . J <„) = -{8x 2 - 18x + 5An + 63)—^ (66) 

Tr(^ 1 [4(2x 2 -6x + 3)c 4+ , +1 ] < l . J < n ) = (2x 2 - 6x + 3){8n(x + A) + 8-^) (67) 
Tr{A- 1 [-27{2x 2 -6x + 3)c l+j }o< t ^<n) = -27(2n + l)(2x 2 - 6x + 3) (68) 

Next we evaluate the traces (|6"T]) - ([6"4"j) . To calculate the first of these (f6T|) We expand the matrix 

A~ l [2{2{i + j) + S)ai +j+1 ] < iJ < n A~ 1 B n 

as 

A' 1 ((4L n + 3I)B n + B n (AL n + 3/)) A" 1 ^ . 
The trace (fBTj) thus can be evaluated as 

Tr((4L„ + 2,I)B n A- l B n A- l )+Tr{A- l B n (AL n + $I)A?B n ) 

=2Tx{A- l B n {AL n + Z^A^Bn) 
=2Tr((4L„ + 3/) {A^Bnf ) 

where we used the symmetry of the matrices in the trace calculation. Using the expression (|23p for 
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A^Bn, we have 



(A^Bn)'' 





1 



<70 OQO n 

(T n _ 2 <7 n _3 + <T n _ 2 C r n 
CT n _l <7 n _ 2 + 0" n _i<7 n 



1 CT r . 



Cn-1 + C„ 



Therefore 



(4L„ + 3J) (A~ l B n y 





1 



3do 

10(72 



3<7 cr n 
7(cj + <?ia n ) 

10((7i + (T2C7n) 



(4n-5)(J n _ 2 (4n - 5)(fJ n _3 + <r n _ 2 (J n ) 

(4n - l)<J n -i (4n - l)(<7 n -2 + ovi-i^Vi) 

1 (4n + 3)a n (4n + 3) (a„_i + a 2 ) 



and thus 



Tr(A" 1 [2(2(i + j) + 3)a i+]+1 ]<^< n A- l B n ) = 2((4n - l)a n _i + (4n + 3)(<t„_i + a 2 )) 



M„ K 2 

4(4n + l)-^ + 2(4n + 3)-^| 



where the last equality is a consequence of the identities (|22p and (|45p . 
Using similar reasoning, we see that the trace term (I62p 

Tr^H&c 2 - 18s + 27(i + j) + SeK+ylcKi^A- 1 ^), 



evaluates to 



3x 2 - 18s + 54n + 36) 



Kn 



The trace term (j63|) is expanded using (|27p . We obtain the following trace identities 



(69) 



(70) 



Ti{A- l [A{2x 2 - 6x + fyci+j^ij^A^-Bn) = 4(2x 2 - 6x + 3)Tr(A; 1 (^ n + A n E r ^)A- 1 B n ) 

= 4(2x 2 - 6x + 3)(Tr{A- 1 E n B n ) + Tr(E^A~ l B n )) 
= 8(2x 2 - 6x + 3) , ft(^'^ 1 B„). 



18 



Since 



E n A n lB n 



oi/2 a 2 /2 
a a x 



the trace reduces as follows: 



Tr(^ 1 [4(2x 2 - 6x + S^+^cKi^A- 1 ^) 



a n /2 a a /2 + ... 
a-n-i a ai + ... 



4(2x 2 - 6x + 3)((2n - l)ai + 2a a n ) 
4(2n - l)(x + 4)(2x 2 - 6x + 3) 
+8(2x 2 -6x + 3)^. 



(71) 



Using the expansion (|30p . the final trace term (|64p simplifies as follows: 



Tr(A- i [-27(2x 2 - 6x + 3)c i+i _i] < J , i <„^- i J B n ) = -54n(2x" - 6x + 3) 
Now we return to equation ([55]) . In this equation, there is a term containing 



(72) 



4Tr(^ 1 S n ) 

We expanded the derivative above as the difference of two traces in ([56]) . We multiplied ([56]) by 
(x — 3)(2x — 3)(4x — 3) and expressed the product 

(x - 3)(2x - 3) (4s - 3)d x Tr(A- 1 B n )) 

as the sum of four terms in ([571) , ([581) , ([59]) , ([60]) minus the sum of four terms in ([6T]) , ([62]) , ([63]) , ([64"]) . 
These eight traces in turn were simplified as ([65]), ([66]), ([67]), ([68]), {SI, ([70]), ([71]), ([72]). Computing 
the sum of ([65]), ([66]), ([67]), ([68]) minus the sum of ([69]), {70]), dTTJ , ([72]) we obtain 



H n (x-3)(2x-3)(4x-3)d x Tr(A- i B n )) = 2(4n + 5)iV n - 2(4n + l)M n - 27K n + 



Ki 



(4s - ll)(2s 2 - 6s + 3)H n - 2(4n + 3)— - 

Now we multiply this through by 2(4n + 3). The left hand side becomes the last term on the left hand 
side of equation (|55p . Substituting and simplifying, we obtain the intermediate identity 



(x - 3) 2 (2x - 3) 2 (4x - 3) 2 d 2 x H n - 
(s - 3) (2s - 3) (4s - 3)(8(n - 3)s 2 - 6(5n - 13)s - 3(9n 2 + 13n + 29))d x H n - 
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2(64nx 4 - 424nx 3 + 2(475n - 6)x 2 - 3(283n - 15)x + 3(91n - Q))H n - 

2(4n + 3)(8nx 2 - 6(5n + l)x - 3(9n 2 + 13n + 17))^+ (73) 

4(4n + l)(4n + 3)M n - 
4(4n + 3)(4n + 5)iV n = 

Recall that equation f|33|) of Lemma [3] expresses d x H n in terms of and -ftT n . Also, Lemma [J] 
gives a linear relationship between ff n , K n , M n and iV n . Thus by using these two equations we can 
eliminate both the d x H n and the M n term in the above formula: 



4(4n + l)(x - l)(x - 3) 2 (2x - 3) 2 (4x - 3) 2 d 2 x H n - 
4(4n + l)(64n(n - l)x 5 - 96(3n 2 - 8n - 2)x 4 + 12(36n 3 + 163n 2 - 65n - 8)x 3 - 
4(459n 3 + 1661n 2 + 949n + 417)x 2 + 3(243n 4 + 2106n 3 + 5192n 2 + 4291n + 1482)x 

-3(243n 4 + 1674n 3 + 3679n 2 + 3140n + 1008))F n + 
8(4n + l)(4n + 3)(16(n + 2)x 3 - 4(17n + 31)x 2 + 6(9n 2 + 48n + 53)x- (74) 

3(18n 2 + 80n + 77))if n - 
32(4n + l)(4n + 3)(4n + 5)(x-l)iV n = 

At this point we have enough information on the interrelationship between H n , K n , N n and M n 
that allows us to derive a number of evaluations of H n (x) for special values of x. These evaluations 
also turn out to be essential for the final step of the proof of Theorem Q] in Section [H 

3 Product form evaluations of H n (x) at special x 

Our calculations for the special values rely on identities we have proved for H n , K n , N n and M n as 
well as a Dodgson-like determinant identity which is useful. This identity is as follows: 

Proposition 1 Suppose the determinants H n , K n , M n , and N n are defined as in t!3\) , \21\) , |3^[ ) ; and 
135\) respectively. Then 

H n ^H n+l = H n N n - H n M n - K 2 n . (75) 

Proof (Proof of Proposition [1]) This identity can be proved using techniques similar to those in [7]. 
We first prove a general determinant expansion result, and then specialize to Hankel determinants to 
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obtain (1751). Consider two matrices 



R"n \j"i,j\o<i,j<.n 
X n +l = [Xi,j]o<i,j<n+l 



where ultimately we will set for all i,j 



Consider the sum 



n+l 


f 






det 


[1 


k=0 







r i,j — 



if i ^ n 



a i+j (x). 



Xkj+i if i = n 



det 



0<ij<n 



if i < C k 
Xi + \ j iii>k 



(76) 



0<i,j<n i 



as a function of the matrix X = X n+ \. It is not difficult to prove that if two adjacent rows of X are 
switched then the sum (|76p changes sign. Since the set of all permutations of the rows is generated 
by flips of adjacent rows, it follows that (|76p is alternating. 

In addition to this (|76p is linear in each row of X. Since (|76p is both alternating and multilinear, 
it is equal to a multiple (depending on R n ) of det(X n+ i). To determine the multiple we set X to the 
matrix: 

1 ... 

r 0i0 r ,i ... r 0j „ 
r 1)0 n ; x ... r 1>n 



X, 



71+1 



r„ i0 r n ,i 



In this case 



det(X n+1 ) = det(R n ) . 

In the sum (j76[) only the term corresponding to i = n + 1 survives and this term itself evaluates to 

det(i? n )det(i? n _i). 

Therefore for all R and X, the sum (|76p evaluates to det(i? n _i) det(X n+ i). When i? and X are Hankel, 
only the last three terms of in the sum (|76p survive. In particular, if we set 



we obtain identity (J75|) . 
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3.1 The evaluation of H n (3) 

We use the identities involving H n , K n , M n , and N n that we already have, along with (175p . At the 
point x = 3 the left hand side of ([75]) evaluates to H n _i(3)H n+ i(3). The right hand side of ([75]) is 
more difficult to evaluate; we have to use the identities from the previous section. At x = 3, equation 
(f33l) simplifies to 

2(4n + 3)AT n = 3(9n 2 + 19n + 14)# n . (77) 

Therefore 

_ 3(9n 2 + 19n+14) 
*»(3) *»(3) (78) 

Now we can use this together with equation (|74p to calculate A n at x = 3 

Ar 3(243n 4 + 1512n 3 + 3605n 2 + 4144n + 1920) TT 

^ (3) = 8(4n + 3)(4n + 5) Fn(3) (79) 

Finally we use the linear relationship of Lemma H] to calculate M n at x = 3. 

3n(243n 3 + 540n 2 + 559n + 58) 
Mn(3) " 8(4n + l)(4n + 3) ^ (3) (80) 



We can now write the right hand side of ()75[) as 

9(3n + 4)(3n + 5)(6n - l)(6n + 1) 2 
4(4n + l)(4n + 3) 2 (4n + 5) n[ ' 

Letting H n = H n (3), we have 



_ 9(3n + 4)(3n + 5)(6n - l)(6n + 1) 2 
n+1 ~ 4(4n + l)(4n + 3) 2 (4n + 5) n { ' 



Writing (|8ip as a recursion for H £ +1 , we find that 



#n+i #i fr 9(3* + 4) (3i + 5) (6i - 1) (6i + 1) 



i=i 

Thus 



tf n # M 4(4i + l)(4i + 3) 2 (4i + 5) 



ff - ff r g Vn Tl ; 9(3^ + 4)(3i + 5)(6i-l)(6z + l) 
n °W /i^y 4(4i + l)(4i + 3) 2 (4i + 5) 

We compute directly that Hq = 1 and i/i = —1. Rearranging 

n ~ ( } jUk JUL 4(4i + l)(4i + 3)2(4i + 5) (82) 
This double product expression for H n can be rewritten in terms of factorials. We record this as a 
corollary. 

Corollary 1 Suppose a^{x) is defined as in ([?j) and H n (x) = det[a,i+j(x)]o<i,j<n- Then 

H /ox_ r . w A (6i-3)l(3i + 2)l(2i-l)l 

ff„(3)-(-l) ll (4i _ 1)!( 4 i + 1)!(3i _ 2)! (83) 
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3.2 The evaluation of #„(§) 

For x = §, the left hand side of identity ([33]) of Lemma [3] vanishes, and we obtain 

3 3(9n 2 + 22n+ll) 3 
*»( 3 ) " H n(~ 2 ) (84) 

Next, we evaluate ([73]) at x = | using this expression for if n . The second derivative terms drops out 
and we obtain 

AT 3 3(243n 4 + 1674n 3 + 4031n 2 + 3934n + 1290) 3 

" 8(4n + 3)(4n + 5) Hn ™ (85) 

Finally, using these expressions for K n {\ ) and iV n (|) in the linear identity ([36]) of Lemma H] at x = |, 
we obtain 

M„(l) = -3"i243n 3 + 702^ + 547,, + ll 8 ) gi 



Substituting these in ([75]) with i/„ = iJ n (|), we obtain 



9(3n + 2)(3n + 4)(6n + l)(6n + 5) 2 

4(4n + l)(4n + 3) 2 (4n + 5) n V 7 

This can be written as a recursion for J? n giving 

tin— 1 

_^ gi jt 1 9(32 + 2)(3z + 4)(6z + l)(6z + 5) 
"~ n-1 "-ffoii 4(4i + l)(4z + 3) 2 (4i + 5) 

Iterating ([55]) . and using the fact that i?o(§) = 1 an d = 2, we obtain the product formula 

„ _on TT W 9 ( 3 * + 2 X 3 * + 4 X 6 * + iKg + 5) 
nU) " Mil 4(4z + l)(4z + 3) 2 (4z + 5) 



The expression (|89[) can equivalently be written in terms of factorials as in (]90p . 

Corollary 2 Suppose a^x) is defined as in ([?]) and H n (x) = det[ai+j(x)]o<ij< n - Then 

H r3 , A («-l)l(6i)l(3i + l) 

^ (2) - 11 2(4i-l)!(4i + l)! (90) 
3.3 The evaluation of H n (\) 

We can mimic these steps for the evaluation of H n at x = | as well. Evaluating the identities (|33p of 
Lemma [3] ([73]) . and ([36]) of Lemma [3] at x = §, we obtain respectively, 

3 3(18n 2 +38n + 19) 3 

*" (s) " 4(4^3) Fn(5) (91) 

Ar „ 3(486n 4 + 3024n 3 + 6724n 2 + 6290™ + 2085) » 

^ = 16(4n + 3)(4n + 5) ^ (92) 

3 _ -3n(243n 3 + 540n 2 + 316n + 31) 3 

n(i) " 8(4n+l)(4n + 3) n(i) 
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Substituting these in ([75]) with H n = iJ n (|), we obtain the recursion 

rr rr _ H3n + l)(3n + 2)(6n + 5)(6n + 7) 2 

n+1 " 4(4n + l)(4n + 3) 2 (4n + 5) " 1 j 

Using the values i?o(f ) = 1 an d -Hi(f) = f > and calculating as in the derivation of (|89p . we obtain 

W jS ,5 4(4 i + l)(4 i + 3rt4i + 5) <94) 
The expression ()94|) can be written in terms of factorials as given in (|95p . 

Corollary 3 Suppose afc(x) is defined as in ([?[) and H n (x) = det[ai + j(x)]o<i t j< n . Then 

3 " (2i-l)!(6i + l)! 
FnW " 11 2(4i-l)!(4t + l)! (95) 

4 A third identity and the differential equation for H n 

We now continue with the proof of Theorem [TJ Using the expression for d\H n in terms of H n , K n and 
N n from (|74p . and the expression for d x H n in terms of and from (|33p 



(a; - 3) (2a; - 3) 2 (4x - 3) 2 ((x - l)(at - Z)d 2 x H n + (2(n + 2) (a: - 3) + 3)^^ - 3n(n + l).ff n ) = 
8(4n + 3)(4n+5)(a;- 1)JV„ 

-4(4n + 3)(-27n 2 - 93n - 75 + (27n 2 + 81n + 60)x + (8n + lO)^ 2 )^ 

+ (-729n 4 - 3564n 3 - 6015n 2 - 4236n - 1080 + (729n 4 + 2916n 3 + 3777n 2 + 1722n + 180)x + (96) 
(432n 3 + 1492n 2 + 1676n + 600)x 2 )H n 

To prove Theorem [TJ we show that the expression on the right hand side of (|96p is identically zero. 
This is the consequence of our third identity. 
Write the right hand side of (j96p in the form 

- (p n N n + q n K n + r n H n ) (97) 

where p n = p n (x),q n = Qn(x), and r n = r n (x) are the negatives of the polynomials multiplying N n , 
K n , and H n on the right hand side of (|96p . respectively. That is, 

p n (x) = -8(4n + 3)(4n + 5)(ar-l) 

q n (x) = 4(4n + 3)(-27n 2 -93n- 75 + (27?i 2 + 81?i + 60).x+(8n + 10)a: 2 ) 
r n [x) = -(-729n 4 - 3564n 3 -6015n 2 - 4236n- 1080+(729n 4 + 2916ri 3 + 3777n 2 (98) 
+ 172271 + 180)a; + (432n 3 + 1492n 2 + 1676n + 600)a; 2 ) 
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Note that p n is linear, q n and r n are quadratic polynomials in x. 

The three matrices in (fT3jh (f2T|) . (f35j) that define /T n , and N n respectively, differ only in their 
last column. Therefore the expression 



p n N n + q n K n + r n H n 
is the determinant of the (n + 1) x (n + 1) matrix whose 
i) first n columns are the columns of the matrix A n in (T[3 



(99) 



ii) whose last column is the linear combination 

PnVn+2 + q n V n +l + r n V n 

where for j > 

aj 
a j+ i 



(100) 



To show that the determinant (|99p is zero, it suffices to show that the last column (jlOOp of this matrix is 
a linear combination of the first n columns vq, v\, . . . , u n _i, i.e. there are weights w n Q, w n i, . . . , w n , n+ 2 
with 



(101) 



w r , 



such that 

n+2 

E Wn.jdi+j = 

i=o 

for i = 0, 1, . . . , n. This is the third identity we need. 



(102) 



After some experimentation, bolstered by checks in Mathematica for up to n = 30, we found that 
the weights are likely given by 

min(-j+n+l, |_f J +l) min(-j+n+2, [_§ J +l) 



Wn,j(x) 



E 

i=0 



i=0 



(103) 



where 
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with 

3 ■ T+HAn + 3)(4n + 5)( 2 ^( J 2 + 2 ^:j+ 1 1 j+ 1 1 )(2n - » + 2)!(4n - 2» + 5)!!(6n + 7)!! 
Un ' 1 {2i + l)(3j + l)(3j + 2)(2i - j - 2n - 3)(2« -j-2n- 2)(2n + l)!(4n + 5)!!(6n - 4i + 7)!! 

a n ,i,j = 3(9j(j + 1) + 2)(12i 2 + 6(3j - 4n - 5)i + 6(n + l)(2n + 3) - j(18n + 19)) 
/3„ i4J - = -3(3j + 2)(24i 3 + 8(12j - 9n - 8)i 2 + 2(36n 2 - 96jn + 64n + 9j(5j - 11) + 23)i 

-2(n + l)(2n + 3)(6n + 1) - f{90n + 89) + 3j(32n 2 + 66n + 33)) 

j ntiJ = 2{2i + 2j - 2n - 3)(i + j - n - l)(36i 2 + 6(9 j - 12n - 5)« + 6n(6n + 5) - 3j(18n + 13) + 2) 

3 • 2 l+3 (2n -2i-2j + 5)(4n + 3)(4n + S)^ 1 ) (^"ij+^i) ( 2n ~ * + 3 ) ! ( 4 ™ - 2i + 5)!!(6n + 7)!! 
Vn,i ~ (3j + l)(3j + 2)(2i -j-2n- 5)(2i -j-2n- A){2i -j-2n- 3)(2n + 2)!(4n + 5)!!(6?i - ii + 9)!! 

<5n,<j = 3(9j(j + 1) + 2)(6i 2 + 3(3j - 4n - 7)t + 3(n + 2)(2n + 3) - j(9n + 14)) 

e 7l) i,j = -3(3j + 2)(12i 3 + (48j - 36n - 50)i 2 + (45j 2 - 3(32n + 49) j + 4(n + l)(9n + 16))i 

-2(n + 2)(2n + 3)(3n + 2) - j 2 (45n + 67) + 3j(n(16n + 49) + 37)) 
e ntiij = 2(2i + 2j -2n-S)(i + j-n- 2)(18n 2 - 9(4i + 3j)n + 33n - 33j + 3i(6i + 9j - 11) + 13) 



The proof of formula ()102f) , and consequently the proof of Theorem [IJ is therefore reduced to a 
somewhat involved system of binomial identities involving summations over multiple indices. Our first 
reaction was that a pen-and-paper proof was infeasible. The explicit form of the conjectured weights 
gave us a slim hope that automated binomial identity provers might work. However, we did not see 
how to do this. So our entire approach to the evaluation of the determinants H n seemed to founder on 
these identities. Our efforts seemed to indicate that these identities may well be beyond the standard 
formulations and computational resources required for automated binomial identity provers, and even 
if such a proof is possible, the insight such a proof would offer would be minimal. 

We were able to find an alternative approach, which avoids the explicit form of the weights, and 
which has the advantage of generalization to other cases. For the proof, all we need to do is to show 
the existence of weights satisfying (|102p and (jlOip . This becomes a manageable problem which relies 
on certain properties of the generating function 

f = f(x,y) = Y,a k (x)y k (104) 

k>0 

of the afc. 
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Lemma 5 Suppose <Xfc(x) is as defined in ^ and f is as in {10J$ - Then 

f( ) = 4x + 2t - 6 

V> (x - 3) (4s - 3)y(2t - 3) - (x - l)(27y - 4) 1 j 



where 



satisfies t 3 y = t — 1. 



The proof can be found in Appendix I. 

The proof of the existence of weights satisfying (|102j) and (llOip takes several steps. Suppose t is 
as in Lemma [5] and define 

r 1 = 2t-3 = -l + 2y + 6y 2 + ••• (106) 

Then 

/<«,») 41 - 3 + ' 



(x - 3) (4a; - 3)yr? - (x - l)(27y - 4) 
The first step is to note that we can find a nontrivial polynomial Qq = Qo(y) of degree n + 1 such that 

vQo = Qi + y 2n+3 ^o (107) 

where Qi = Qi(y) is a polynomial of degree n + 1, and = ^o(y) is a power series in y. This is 
because there are n + 2 coefficients in Qo and only n + 1 linear equations these coefficients must satisfy. 
In the next step, put 

Q 2 = [x - 3) (4s - 3)yQ 1 - (x - l)(27y - 4)Q (108) 

Then Q2 = Q2(x,y) is a polynomial in x and y of y-degree n + 2. It is important to note that all 
three of the polynomials Qq, Qi, Q2 so defined are nontrivial. We claim that the coefficients of Q2 are 
the weights we want. In other words, the coefficients of the terms y n+2 through y 2n+2 in /Q2 vanish. 
Our aim now is to use the special form of the generating function / in (|105|) to show this. 
Writing f)105j) in the form 

f(x, y)((x - 3)(4x - 3)yrj - (x - l)(27y - 4)) = 4x - 3 + 77 

and multiplying through by Qo, we get 

f(x, y)((x - 3) (4c - 3)y(Qi + y 2n+3 ^o) - (x - l)(27y - 4)Q ) = (4x - 3)Q + Qi + 2/ 2n+3 ^o 

or 

fix, y) (Q 2 + (x - 3)(4x - 3)y 2n+4 V ) = (4x - 3)Q + Qi + y 2n+3 *o 
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and therefore 

f(x, y)Q 2 = (4s - 3)Q + Qi + V 2n+3 (*o + /(*, y)(x - 3) (4c - 3)j/*„) 
which means that 

/Q 2 = (4x - 3)Q + Qi + y 2n+3 ^i (109) 

where 

(4s - 3)Q + Q\ 

is a polynomial in x and y of y-degree n + 1 and = ^f(y) is a power series in y. 
This last statement (|109|) is equivalent to the statement that 

n+2 

J2 C n+2-j(Q2)a i+j = (110) 

i=o 

for i = 0, 1, . . . , n, where by Cfe(^) we denote the coefficient of the term y k in a power series Thus 
([102]) holds with 

W nj ' = C n+2 -j(Q 2 ) 

for j = 0,1,..., n + 2. 

In particular, taking determinants, we have that 

C {Q 2 )N n + Ci(Q 2 )# n + C 2 (Q 2 )H n = (111) 

This identity is not trivial. Otherwise, we would have a nontrivial linear relationship among the n + 1 
columns Vq, v\, . . . , v n of H n . But we have already shown that the special values H n (3) , H n (^) , H n (j) 
evaluated in ([83]) . ([90]) . ([95]) are nonzero, and therefore H n does not identically vanish. 

Now we need to verify the three special values in (jlOip . Since equation (jllOp is valid when 
multiplied through by an arbitrary constant, in terms of the coefficients Ck(Q 2 ), we need to show that 
for some nonzero constant a, 

Co(Q 2 ) = aPn 

Ci(Q 2 ) = aq n (112) 
C 2 (Q 2 ) = ar n 

where p n ,Qn, r n are as defined in ([05]) . 
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First we rewrite (jllip in terms of the three coefficients Cq(Qo),Ci(Qo),C2(Qo) which are pure 
constants, independent on x and y. First, comparing coefficients in (|107j) and (|108p . we obtain the 
five identities in (11130 . 

Co(Qi) = -Co(Qq) 

C x (Qi) = 2C (Q )-Ci(Qo) 

C (Q 2 ) = 4(x - l)Co(Qo) (113) 

C X (Q 2 ) = (x-3)(4x-3)C (Qi)-27(x-l)C (Qo) + 4(x-l)Ci(Qo) 

C 2 (Q 2 ) = (a;-3)(4x-3)Ci(Q 1 )-27(x-l))Ci(Qo) + 4(x-l)C 2 (Qo) 



Therefore 



C (Q 2 ) = 4(x - l)Co(Qo) 

Ci(Q 2 ) = -2(22; 2 + 62;-9)Co(Qo) + 4(x-l)Ci(Qo) (114) 
C 2 (Q 2 ) = 2(a;-3)(4a;-3)Co(Qo)-2(2x 2 +6x-9)Ci(go) + 4(a;-l)C 2 (Qo) 

and (jllip becomes 

4(x-l)C (Q )N n + 

(-2(2x 2 + 6x - 9)C (Qo) + 4(x - l)Ci(Qo)) + (115) 
(2(x - 3)(4x - 3)C (Qo) - 2(2x 2 + 6x - 9)d(Q ) + 4(x - 1)C 2 (Q )) H n = 

We have explicit linear relationships between the three determinants N n , K n , H n at the points x = 
3, |, |. For x = 3, we use the expressions in terms of H n (3) for iV n (3) from ([79^1 and for if„(3) from 
(|79p in (|115p to obtain the linear equation 

8(3 + 4n)(5 + 4n)C 2 (Q ) + 

6(5 + 4n)(l + 2n + 18n 2 )Ci(Q ) + (116) 
3(30 + 67n + 338n 2 + 540n 3 + 243n 4 )C (<2o) = 

For x = |, we use the expressions in terms of ff n (|) f° r -^n(f) from ([85]) and for if n (§) from ([5"4"|) in 
(fTT5j) to obtain 

8(3 + 4n)(5 + 4n)C 2 (Q ) + 

12(5 + 4n)(2 + 10n + 9n 2 )C 1 (Q ) + (117) 
3(120 + 778n + 1445n 2 + 1026n 3 + 243n 4 )C (Qo) = 
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Finally for x = §, we use (|92]) and (JH]) in (fTToT) to obtain 

8(3 + 4n)(5 + 4n)C 2 (Q ) + 

12(5 + 4n)(-4 + n + 9n 2 )Ci(Q ) + (118) 
3(-240 - 446n + 95n 2 + 540n 3 + 243n 4 )C (Qo) = 

Using (any two of) these equations we obtain the parametric solutions 

C 2 (Qo) = 3(50n + 343n 2 + 540n 3 + 243n 4 )a 

d(Q ) = -12(l + n)(30 + 67n + 36n 2 )a (119) 
C (Qo) = 8(15 + 32n + 16n 2 )a 

Going back to the system (|114jl . we compute that 

C (Q 2 ) = 32(4n + 3)(4n + 5)(x-l)a 

Ci{Q 2 ) = -16(4n + 3)(-75 - 93n-27n 2 + 60x + 81nx + 27n 2 a; + 10x 2 + 8na; 2 )a (120) 
C 2 (Q 2 ) = 4(-1080 -4236n-6015n 2 - 3564n 3 - 729n 4 + 180x + 1722nx + 3777n 2 x + 
2916n 3 x + 729n 4 x + 600x 2 + 1676nx 2 + 1492n 2 x 2 + 432n 3 x 2 )a 

Taking a = —\, these are exactly the weights w n ^ n+ 2 = p n , Wn,n+i = Qn, Wn,n = r n as claimed in 
(llOip . This finishes the proof of Theorem [TJ • 

4.1 The differential equation solution for the (3, l)-case 

In the previous section we derived the differential equation (fT4"j) satisfied by 



a a± a 2 ■■■ a n 
a\ a 2 03... a n+ i 

H n (x) = det 



«2r, 



where 



a k = a k' 1] ( x ) = J2 



m=0 



3k + 1 — m 
k — m , 



We obtain a power series solution to (|14p of the form 



y(x) = 5>(x - 3)* 

i=0 
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where the satisfy 



(3n + i + 3)(n - i) 
Oj+l = — - w„, — — Oi 



(fc + l)(2A: + 3) 



with 6q = 1. This allows us to prove that 



H n (x) = C n J2 



n\{3n + i + 2)\2 i {x-3) i 



(121) 



^(3n + 2)!(n-i)!(2i + l)! 

where C n is the constant of integration that depends only on n but not on x. We have already 
calculated C n = H n (3) in (|83|) of Corollary [TJ We record the final determinant evaluation result of 
the (3, l)-case as a Theorem: 



Theorem 2 Suppose ak{x) is defined as in ftTty and H n (x) = det[ai + j (x)]o<i,j<n- Then 

H (r) _ ( 1 \n TT ( 6 * ~ 3 )K^ + 2)!(2i - 1)! " n!(3n + i + 2)!2' (x - 3)' 
n t ii (4* - l)!(4i + l)!(3i - 2)! ^ (3n + 2)!(n- *)!(2» + l)!" 

In particular, evaluating (|122p at a; = and at x = 1 we obtain 
Corollary 4 



(122) 



det 



dct 



3(i + j) + l 
* + .7 

3(i + j)+2' 
i + j 



f n »n (6i ' 3)!(3 ' + 2)!(2i ' 1)! v 

1 J 11 (4i- 1)1(4* + l)!(3t- 2)1 

(6i - 3)!(3i + 2)!(2i - 1)! <A n!(3n + i + 2)!(-4) 



n!(3n + si + 2)!(-6) J 



0<i,j<n 



L (4i - l)!(4i + l)!(3i - 2)! ^ (3n + 2)!(n - i)!(2i + 1)! 



■E 



Note that the second evaluation (|124p is an alternate form of the product formula ([3]) . 



(123) 
(124) 



4.2 A recursion for H^'^( 



x 



We now remark on a property of the polynomials H n (x) for the (3, l)-case, which does not enter into 
the evaluation of the determinant, but which nevertheless is of interest. 

Theorem 3 With Hq = 1, H\{x) = —2x + 5, the H n (x) satisfy a three-term polynomial recursion 



p n {x)H n (x) + q n {x)H n ^ 1 (x) + H n _ 2 (x) = 



(125) 



for n > 2 where 

Pn(x) = 
Qn(x) = 



4(4n-3)! 2 (4n- l)! 2 



9(3n - 2) 2 (3n - l) 2 (2n - l)! 2 (6n - 5)! 2 
4(n - l)(4n - 5)!(4n - 1)! 



(x-1) 2 

(8(4n-3)(4n + l)x 3 



3(3n - 2)(3n - l)(4n + l)(2n - l)!(6n - 5)! 

-36(4n - 3)(4n + l)x 2 + 6(126n 2 - 63n - 23)x - 3(108n 2 - 54n - 19)) 
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Proof The theorem follows from the explicit form of the H n (x) given in Theorem [2j We omit the 
details. • 



The following property of the H n for the (3, l)-case is a consequence of Theorem [3) 
Corollary 5 The roots of H n {x) are real and interlace those of H n+ \{x) for n > 0. 

Proof The polynomials H n (x) are clearly not orthogonal, so the proof of the interlacing phenomenon 
does not follow from classical results, nor directly from the results of [8]. Nevertheless an elementary 
induction argument works, which we also omit. • 

It is also curious that the recursion (|125p seems to define a sequence of rational functions with 
nontrivial denominators. We are unaware of a general theory of such recurrences which actually define 
Sturm sequences of polynomials. 

Another consequence of the recursion in Theorem [3] is a completely different derivation of the 
product formula ([3]) for the (3, 2)-case. Letting x = 1 in (|125p . and noting p n (l) vanishes, we obtain 
the recursion 

#(3,2) = - 1 #(3,2) 

3(3n + l)(3n + 2)(2n + l)!(6n + 1)! (3 ,2) 



4n(4n + l)(4n- l)!(4n + 3)! 

from which the product evaluation 



H n-1 



(3i2) = fr 3(3i + l)(3z + 2)(2z + l)!(6z + l)! 
f = \ 4t(4i + 1)(4« - l)!(4i + 3)! 



immediately follows. 

5 Transformation rules 

(2 1) 

Before we proceed with the evaluation of the determinant H\ ' (x), we would like to collect the 
information we already have about various trace calculations to automate this process to some extent. 
In the course of the calculations involved in the (3, l)-case, we have had to evaluate traces such as 

Tr ( j4 n 1 b( a; ) a l +j+i]o<i,j<n), Ti(A~ l [q(x){i + j)a i+j+1 ]o<i,j< n ), Ti(A^ 1 [r(x)c i+j } <ij< n ) 

for polynomial coefficients p(x),q(x),r(x). These traces are reduced to calculating traces of the form 

Tt(A^ 1 [a i+j+1 ]), Tr(A- 1 [(i+i)a i+ ,- +1 ]), Tr (A^fcj+j]) 
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etc., and then use linearity of the trace. Furthermore, the calculations do not use the explicit form 
of the polynomials a^. We can go one step further: when we apply the operator Tr^" 1 *) to various 
terms of an identity linearly, we can then multiply through the resulting transformed identity by H n 
to obtain a relationship between various determinants related to H n , as we did in a number of cases 
in the computation of the (3, l)-case. In Table [1] we summarize the effect of this transformation on 
various terms that may appear in an identity. 











na„_i — > 


c„-i -> 


a n ^ (n + l)H n 


na n — >■ n(n + l)H n 


c n -> (2n + l)a i? n 


a n+l ~^ K n 


na n+ i -> 2nK n 


c n+ i ->■ 2a -fC n + 2naiH n 


a n+2 ->■ M n + N n (n > 0) 


na n+2 -> 2(n - l)M n + 2nN n (n > 0) 





Table 1: Some linear transformation rules. 

In TableHJ a n = a n (x),n = 0, 1, ... is a sequence of polynomials, c n are the convolution polynomials 
defined by 

n 

C n = ^2 a i a n-i (126) 
8=0 

and H n , K n , M n , N n are the following (n + 1) x (n + 1) determinants defined by the a^'s: 



H„ = det 



ao a\ a,2 
ai a 2 a 3 



a n +i 

&2n 



K„ = det 



(l0 CL\ 
O'n—l 0>n 



On-1 On+1 

d n CLn+2 

«2n-2 «2n 

«2n-l 02n+l 



Mr, = det 



ao ai 

a\ CL2 

O'n—l 

0"n Q"n+1 



On-2 On 

«n-l On+2 On+1 

fl2n-3 «2n «2n-l 

0-2n-2 0-2n+l &2n 



N n = det 



a\ CL2 



«n-2 On-1 On+2 



On+3 



fl2n-3 «2n-2 fl2n+l 
«2n-2 «2n-l ^2n+2 
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From these rules we obtain a few known applications that are not connected with the determinants 
we are evaluating. We give a few simple examples. 

Example 1 Recall that a polynomial sequence {a,k}k>o is sa id to be an Appell set if d x a n = na n -\. 
If B n , 7i n , T n , P n are the Bernoulli, Hermite, Chebyshev, and Legendre polynomials, respectively, 
then for any complex a / 0, 

{B n (x)}, {(2a)-«H n (ax)}, {(x 2 -a 2 )%T n ( X )}, {(x 2 - a 2 )% P n ( % )}, 

Vi — or \l x L — a z 

are all well-known examples of Appell sets [HE]- Starting with the defining identity 

d x a n = na n _i 

and applying the transformations d x a n — ► d x H n and na n -\ — * from Tabled! we see that d x H n = 0, 
and therefore the Hankel determinants of Appell polynomials are independent of x. 

Example 2 For a sequence of polynomials a n satisfying 

xd x a n = na n + cna n -i (127) 

for some constant c, W6 make th.e replacements d x a n —> d x H n , na n — > n(n + l)il n and na n „i — > 
from Table Q] to find that the corresponding Hankel determinant H n satisfies 

xd x H n = n(n + l)H n 

This is integrable with 

Hfi{x^j — c n x 

For c = —1 in (11271) we get the Laguerre polynomials L n with Lq = 1,L\ = 1 — x. After simple row 
operations, c n = H n (l) has a known product evaluation as 



Cn = det 



1 



2- n2 (-i)^n 1 



0<„<n ^ ' ^■^■■■(2i-l) 



n—i 

! 



.(i + i)! 

For c = 1 in (|127p we get the so-called derangement polynomials 

n I 

. 77 ! 

i=0 

for which the constant of integration is 

n 

c n = Y[n 2 

8=1 

as found in |10|. 
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Example 3 Legendre polynomials satisfy 

(x 2 -l)d x P n = nxP n - nP n -i 

with Po = 1, Pi = x and therefore 

(x 2 — l)d x H n = xn(n + l)H n 

and 

H n = c n - (x 2 - 1) 2 . 



i/ n can be evaluated at x = 0, giving 



,2 



Example 4 In general, for any family of polynomials {ak}k>o such that 

satisfy a linear identity with fixed polynomial coefficients (where c n is the convolution defined in (|126p ). 
the corresponding Hankel determinant is the solution to a simple separable first-order linear equation. 

Now we return to the determinants H n ^' a ^ and consider the (2, l)-case. 



\x m 



6 The (2, l)-case 

Now we let 

(2,1)/ v A /2fc + 1 - m N 
ak = a * {X)= ^A k-m , 

m=0 \ / 

k 

Ck = ^2aia k -i (128) 
i=0 

H n = H n 2 ^(x) = det[a i+: ,-(x)]o<i,j<r i . 

(2 1) (2 1) r— — 

Theorem 4 Let the polynomials df. = a k ' (x) and H n = H n ' (x) are as defined in 11128]) . Then H n 
satisfies the differential equation 

x(x - 2)d 2 x y + (2x - l)d x y - n(n + l)y = (129) 

The proof of Theorem [5] mimics the steps of the (3, l)-case. The analogues of the first two identities 
in Lemma [T]and Lemma [2] of the (3, l)-case are now as follows: 
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Lemma 6 Suppose and are as defined in t!28\) . Then 

2x(x - 2)d x a n - (n + l)a n+1 + (An + 2(x + l))a n - (x - l)c„ + 4(x - l)c n _i = (130) 

Lemma 7 Suppose ak and Ck are as defined in M28\) . Then 

(nx + 2(x + l))a n+2 - (2x(x + 2)n + 2(2x 2 + 3x + 4))a n+ i 

+4x 2 (2n + 3)a n + (x - l)(x - 2)c n +i - 4(x - l)(x - 2)c n = (131) 

The proofs are given in Appendix II. 

Applying the transformations in Table [Tj to the identities in these two Lemmas, we immediately 
obtain the pair of identities 

2x(x - 2)d x H n - (2n + l)K n + 4n(n + l)H n + 2(x + l)(n + 1)#„ - (x - l)(2n + l)H n = 
x(2(n - l)M n + 2n7V n ) + 2(x + l)(M n + iV n ) - 2x(x + 2)(2nK n ) 
+2(2x 2 + 3x + 4)K n + 8x 2 n(n + l)H n + 12x 2 (n + l)JET n 

+(x - l)(x - 2)(2K n + 2n(x + 3)H n ) - 4(x - l)(x - 2)(2n + l)H n = 

These can be rewritten as 

2x(x - 2)d x H n + (x + 3 + 8n + 4n 2 )F n - (2n + l)K n = 0, 

(2nx 3 + (8n 2 + 12n + 8)x 2 + (12n + 10)x -An- 8)H n (132) 
-((2 + 4n)x 2 + (12 + 8n)x + A)K n + 2(nx + l)M n + 2((n + l)x + l)JV n = 

First we use identity (|25p to rewrite K n as a trace in (|132|) . Differentiating, we follow through 
calculations of traces similar to the derivation of d x H n of the (3, l)-case to obtain 

2(x - 2) 2 x 2 (l + nx)d 2 x H n 

+2(x - 2)x(3 + lOn + 4n 2 + 3x + 9nx + 12n 2 x + 4n 3 x + 4nx 2 + 2n 2 x 2 )d x H n 
+ (10 + 49n + 78n 2 + 44n 3 + 8n 4 + (4 + 34n + 86n 2 + 98n 3 + 48n 4 + 8n 5 )x 

+ (1 + 8n + 14n 2 + 8n 3 )x 2 + (2n - 2n 3 )x 3 )H n (133) 
-(2n + 1)(1 + 2n + 2nx + 2n 2 x)N n = 

Combining (|133p and the expression (|132p for d x H n , we obtain 

2(x - 2)x(nx + 1) (x(x - 2)d 2 x y + (2x - l)d x y - n(n + l)yj = (134) 
(2n + 1)(1 + 2n + 2nx + 2n 2 x) (7V n - (4 + 2n + x)K n + (2 + 5n + 2n 2 + 3x + 2nx)F n ) 
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6.1 Product form evaluations at special x for the (2, l)-case 

At this point, we can evaluate H n (x) at special points x easily. At x = 0, we get from (|132p and (1133ft 

3 + 8n + 4ra 2 
2n + 1 



10 + 49n + 78n 2 + 44n 3 + 8n 4 
M n = -(3n + 2n 2 )H, 



N " = (2^TTP H " (135> 



Using (fT5l) and that #0 = = 1 f° r £ = 0, we obtain the recursion for J* n as 

n — 1 



(136) 



and therefore 



as is well known. At x = 2 



H n (0) = 1 (137) 



5 + 8n + 4n 2 
* n = 2n+l Hn 

22 + 33n + 20n 2 + 4n 3 

^ " 2^TI F " (138) 

n(7 + 8n + 4n 2 ) 

M n = — -— ti n 



and H = I, Hi = -3. From ([75 



_ (2n-l)(2n + 3) rr2 

^n+l^n-l - (2n + 1) 2 

and therefore 

ff n (2) = (-l) n (2n + l) (139) 

After the evaluations at special points, we now return to the proof of the differential equation for 
the 2n + 1 case. By (|134j) . it suffices to prove the identity 

N n - (4 + 2n + x)K n + (2 + 5n + 2n 2 + 3x + 2nx)H n = (140) 

This is the determinantal form of the third identity for the (2, l)-case. The left hand side of (|140p 
is a determinant as in (|99|) and (|100p of the (3, l)-case. The problem is again to show the existence 
of weights tt> n ,0j • • • ,Wn,n+2 satisfying the third identity (jlOip and (|102p where in this case the 
requirements are 

W n ,n+2 = Pn{x) = 1 

w n , n+1 = q n (x) = -(4 + 2n + x) (141) 
Wn,n = r n (x) = 2 + 5n + 2n 2 + 3x + 2nx 
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By experimentation, we found that these weights are explicitly given by 

^Hi-j;:) + (:-j;3 + «"^r) + (:-^:))^} 

for i = 0, 1, . . . , n— 1, and Wn^+ijWn.n+a as in (|141h . The resulting identities are simpler than the 
(3, l)-case, and in fact a combination of paper-and-pencil and automatic identity provers can be used 
to prove that these are correct. However we can do better, and we use this case as another exercise 
for the application of the method of the existence of the weights that we used in the (3, l)-case. 
We need a lemma similar to Lemma [5l 

Lemma 8 Suppose ak(x) is as defined in (0|) and f is as in 1 104\ )- Then 

K-ti- f.x-sJ+i-^ (143> 

where 

satisfies t 2 y = t — 1. 



The proof of Lemma can be found in Appendix I. 

Let t be the power series in y given in Lemma [H Then there is a nontrivial polynomial Qo = Qo(y) 
of degree n + 1 such that 

tQo = Qi + y 2n+3 *o (144) 

where Qi = Qi(y) is a polynomial of degree n + 1, and = ^o(y) is a power series in y; i.e. the 
coefficients of y k in tQo vanish for n + 2 < k < 2n + 2. Such a nontrivial Qo exists because there are 
n + 2 coefficients and only n + 1 linear equations to satisfy. In the next step, put 

Q 2 = (x - 2)yQ 1 + (1 - 2xy)Q (145) 

Then Qi = Q2(x,y) is a polynomial in x and y of y-degree n + 2. All three polynomials Qo,Qi,Q2 
are nontrivial. We claim that the coefficients of Q2 are the weights we want. In other words, the 
coefficients of the terms y n+2 through 

y 2n+2 in 

/Q2 vanish. Writing (11431) in the form 
f{x,y)((x-2)yt + l-2xy) = t 
and multiplying through by Qo, we get 

f(x, y)((x- 2)y(Q 1 + y 2 " +3 * ) + (1 - 2xy)Q ) = Qi + y 2n+3 *o 
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or 

fix, y)Q 2 = Qi + y 2n+3 *o - /(*, y)(x - 2) yy 2n+3 * 

and therefore 

f(x, y)Q 2 = Qi + y 2n+3 (*o - f(x, y)(x - 2)y^ ) 

which means that 

fQ 2 =Q 1 + y^+H 1 (146) 
where = ^(y) is a power series in y. This last statement (1146P is equivalent to 

n+2 

E C "+2-,(Q 2 )a l+i =0 (147) 

3=0 

for i = 0, 1, . . . , n, where by Ck{^) we denote the coefficient of the term y k in a power series Thus 
(fT02l) holds with 

Wn,j = C n+2 ^j(Q 2 ) 

for j = 0,1,..., n + 2. 
Therefore 

C (Q 2 )iVn + C x (g 2 )K n + C 2 (Q 2 )H n = (148) 

This identity is not trivial, for otherwise we would have a nontrivial linear relationship among the 
n + 1 columns Vo, vi, ■ ■ ■ , v n of H n . But we have already shown by the evaluations in (|137p . (| 139j) that 
H n does not identically vanish. 

Now we need verify the three special values in (|14ip . i.e. for some nonzero constant a, 

Co(Q2) = ap n 

Ci(Q 2 ) = aq n (149) 
C 2 (Q 2 ) = ar n 

where p n , q n ,r n are as defined in (|14ip . 

Rewrite (I148p in terms of Co(Qo) > Ci(Qo) > C 2 {Qq) which are pure constants, independent on x and 
y. Comparing coefficients in (|144p and (|145p . we obtain 

Co(Qi) = Co(Qo) 

Ci(Qi) = C 1 (Qo)+C (Q ) 

C (Q 2 ) = Co(Qo) (150) 

Ci(Q 2 ) = (x-2)Co(Q 1 )-2xC (Qo)+C 1 (Qo) 

C 2 {Q 2 ) = {x-2)C 1 (Q 1 )-2xC 1 (Q )+C 2 (Q ) 
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Therefore 



C (Q 2 ) = C (Qo) 

Ci(Q 2 ) = C 1 (Q )-(x + 2)C (Q ) (151) 
C 2 (Q 2 ) = C 2 (Q )-(x + 2)C 1 (Q ) + (x-2)C (Q ) 



and (j!48j) becomes 



C (Qa)N n + 

(C 1 (Qo)-(x + 2)C (Qo))^n+ (152) 
(C 2 (Qo) - (* + 2)Ci(Q ) + (x - 2)C (Qo)) # n = 

Using the expressions for if n and iV n in terms of H n at x = and x = 2 from f)135 j) and p38|) 

(10 + 49n + 78n 2 + 44n 3 + 8n 4 )C (Q ) 

+(3 + 8n + 4n 2 )(2n + \){C x {Qq) - 2C (Qo)) 

+(2n + 1) 2 (C 2 (Q ) - 2d(Q ) - 2C (Q )) = 
(22 + 33n + 20n 2 + 4n 3 )C (Q ) 

+(5 + 8n + 4n 2 )(Ci(Q ) - 4C (Q )) 

+(2n + l)(C 2 (Qo)-4Ci(Q )) = 



The parametric solutions are 



C 2 (Qo) = n(2n + l)a 

Ci(Qo) = -2(n + l)a (153) 
Co(Qo) = a 



Going back to the system (|151j) 



Co(Q 2 ) = a 

Ci(Q 2 ) = -(4 + 2n + x)a (154) 
C 2 (Q 2 ) = (2 + 5n + 2n 2 + 3x + 2nx)a 

With a = 1, these are exactly the weights u> n ,n+ 2 = Pn, Wn,n+i = ^n, Wn,n = r n as claimed in (j 1 4 1 [) . 
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6.2 The differential equation solution for the (2, l)-case 



We have the differential equation (1129|) satisfied by the Hankel determinants H n for the (2, l)-case. 
We obtain a power series solution to ()129[) as 

" (n + Q!y(x-2)» 



where he constant of integration C n = H n (2) is given by (|139p . Therefore 

Theorem 5 Suppose dk{x) is defined as in |2j) and H n (x) = H n ' (x) = det[ai+j(x)]o<ij< n - Then 



H n (x) = {-l) n {2n + l)Y J 



(n + t)!2*(a; - 2)* 
- (n-*)!(2i + l)l 



(155) 



Evaluating (|155p at x = and at x = 1 we obtain 
Corollary 6 



det 



det 



2(i+i) + r 

'2(i+j) + 2 x 



-l)»(2n + l)£ 



0<i j'<n 



(ra + i)!(-4) i 
-(n-i)!(2i + l)! 



(-lH2n+l)£; 



(n + i)!(-2) i 



(156) 
(157) 



.r;(n-t)!(2i + l)! 

- 0<ij<n «=o v ' v ' 

Interesting point about Corollary [6] is that the right hand side of (|156f) is a complicated way of writing 
1, while the right hand side of (|157p evaluates to the simple expression 



n(n + l) 
(-1)* 



(158) 



Alternate derivations of the evaluation (|158p for the determinant in (|157|) were communicated to us 
by Ira Gessel and Christian Krattenthaler. A generalization of the determinant in (|157p can be found 
in Corollary [8] below. 

(2 1) 

It is also interesting that the determinants H n ' (x) are orthogonal polynomials: 

(2 1) 

Corollary 7 The (n+ 1) x (n + 1) Hankel determinants H n (x) = H n ' (x) /orm an orthogonal family 
of polynomials. 

Proof We have 

H (x) = l, H 1 (x) = l-2x 

and for n > 2, 

H n (x) = 2(1 - x)fl„-i(s) - fl„-2(a:) (159) 

which can be verified directly from the explicit formula in (|155p or using the differential equation 
(I129|) , Orthogonality now follows from Favard's theorem. • 
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In this case the generating function of the H n {x) is 

; ~ 2 = E H nW 

I - 2(1 -x)y + y 2 ^ Q 
Comparing with the generating function 

1 



E U n {x)y n 



l-2xy + y 2 n > Q 
of the Chebyshev polynomials of the second kind U n (x), we have 

H n {x) = U n (l - x) - ET n _i(l - x) . (160) 

In fact we can prove more. From the identity 

m=0 \ / m=u \ / 

and the differential equation fll29|) for Hn ' (x), we obtain both the evaluation of the Hankel deter- 
minant of the polynomials 

M*)=£(f + V, (162) 

m=0 \ / 

(2 1) 

and an alternate expression for Hfi"(x) itself. 

Corollary 8 Suppose the polynomials ak(x) are as defined in /tl6$fy and H n (x) = det[oi+ 3 -(a;)]o<t,j<fi' 
Then H n (x) satisfies Jacobi's differential equation 

(x 2 - l)d 2 x y + (2x + l)d x y - n{n + l)y = (163) 

and therefore in terms of the Jacobi polynomials 

(Vi,n2 , l _u 

H n (x) = {-AT^LP^' »>(*). 

sn(n\) 2 D (±,-|), 



7 Remarks 



When we started to look at Hankel determinants of A n = \<ii+j\o<ij<n with entries = ( 3fc ^ 1 ), it 
looked like an almost product could account for the combination of small and large prime factors that 
were showing up in the data. At the time, we also thought that whatever was moving the (3, l)-case 
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off the pure product formulas of the (3, 0) and (3, 2)-cases would also explain the (3, 3)-case, and the 
(3, 4)-case, and so on. It might also get at the (4, l)-case, and all the rest of the (J3, a)-cases. 

This section is a brief discussion of what did and did not work out. We started with the observation 
that if you define ak(x) by (J7|) then 



3k + 2 X 
. k , 



so that 

fljw) = fljw)(i) 

(3 1) 

has a product form evaluation. What we wanted was the value of Hn ' (x) not at x = 1 but at x = 0. 
Experimenting with the expansion ofHr '(x) at various points, we found that at x = 1 and at x = 3 
the Taylor coefficients of H n ' (x) had small prime factors. We called x = 1 and x = 3 round points. 

We also found several Hankel determinants that had similar behavior but we had trouble finding a 
consistent explanation for this phenomenon. In others, such as H n 3 ' ^ , the coefficients were somewhat 
round but not round enough to allow us to find a simple power series expansion. 

When we tried to determine the common elements of these different determinants differential 
equations came to mind because hypergeometric series such as ([8]) clearly have a differential equation. 
In retrospect, the explanation of why x = 1 and x = 3 are a good places to expand H n ' is partially 
explained by examining the form of its differential equation. For instance (|14p can be written in the 
form 

(x - 3) 2 d 2 x y + 2(n + 2)(x - 3)d x y - 3n{n + l)y = -2(x - 3)d 2 x y - 3d x y. (164) 

This form is significant because the operator on the left hand side of (I164j) takes (x — 3) k to (x — 3) k 
times a polynomial in n and k and the operator on the right takes (x — 3) k to (x — 3) fc_1 times a 
polynomial in n and k. This means that the differential equation defines a two-term recursion on 
successive coefficients of the expansion of Hn ,l \x) around (x — 3). So the power series is easily shown 
to be hypergeometric. 

Differential equations can also explain the "near" round behavior seen in some cases, such as the 
Hn' ■ However, the existence of a differential equation would not provide any input on why the 
constant of integration such as Hn'°\3) was round. Consequently it was a nice bonus when, during 
the proof of the differential equation, we realized that we could prove product formulas for these kinds 
of determinants (e.g., (JS3} in Corollary [TJ (|SQJ) in Corollary El (ESJ) in Corollary El ([Wjl . (fTHHI) . and 
(jTTTI) in Section EJ). 

Once differential equations entered the picture, it was easy to guess the differential equations for 
many of these Hankel determinants and obtain strong experimental evidence for their correctness, but 
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2(3x - l) 2 (x - 3) 3 x(4(n + 2)(2n + l)x 2 + (8n 2 + 20n + ll)z - 1) d*y 
+ (x - 3) 2 (3x - 1) (12(n + 2)(2n + l)(8n + 27)x 4 - 3(128n 3 + 568n 2 + 724n + 161)x 3 

- (576n 3 + 3016n 2 + 4756n + 2269)x 2 + (72n 2 + 252n + 319)x - 15) d\y 
+ 3(z - 3) (12(8n 4 + 118n 3 + 427n 2 + 533n + 174)x 5 

- (736n 4 + 6704 + 19628n 2 + 21289n + 5814)x 4 
+ (800n 4 + 2944n 3 + 564n 2 - 7580n - 7078)x 3 

+ 2(816n 4 + 6744n 3 + 19358n 2 + 23069n + 9809)x 2 

- 6(108n 3 + 540n 2 + 972n + 679)x + 15(9n + 20))d 2 y 

- 3 (12(16n 5 + 62n 4 - 7n 3 - 293n 2 - 378n - 120)x 5 

- (960n 5 + 3104n 4 - 3080n 3 - 20993n 2 - 23419n - 6450)x 4 

+ 4(144n 5 - 4n 4 - 2300n 3 - 4657n 2 - 2040n + 838)x 3 

+ (1728n 5 + 6624n 4 - 5832n 3 - 556990n 2 - 80626n - 36004)x 2 

+ 12(318n 4 + 1995n 3 + 4670n 2 + 4928n + 2058)x - (783n 2 + 3105n + 3102))4y 

- 3n(n + 1) (12(12n 4 + 68n 3 + 137n 2 + 113n + 30)x 4 

- (864n 4 + 4488n 3 + 8158n 2 + 5887n + 1222)x 3 
+ (720n 4 + 2304n 3 + 280n 2 - 4547n - 3388)x 2 

+ 3(576n 4 + 3816n 3 + 9182n 2 + 9533n + 3666)x - 3(120n 2 + 507n + 538))y = 

Figure 1: Differential equation for Hn ' {%)■ 

we did not find a general pattern, even for the families . It is doubtful that such a pattern 

exists. For example, the (unproven) differential equation for the (3, 2)-case is fourth order and quite 
mysteriously complex (see Figure [1]) . This differential equation has been tested for Hn ' (x) where 
n = 1,2, ... , 75. We did not find a differential equation for (3, 3) or for any (3, a) with a > 3. 

(3 1) 

Now the actual proof of the differential equation for a determinant like H n ' (x) is an elusive 
creature. The straightforward evaluation of the derivative yields a sum of n + 1 determinants, each 
one of which is badly behaved at x = 3. That was a disturbing development, and we still had to 
differentiate each of these determinants themselves to get at the second order differential equation. 

In order to approach the problem systematically, we wrote the derivative of a Hankel determinant 
H n = det[aj + j(x)]o<i,j<n as 

d x H n = TriA^dxA^Hn 
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where 

d x A n = d x A n (x) = [d x a i+j (x)}o<i,j< n - 
The obvious thing to try at this point is to look for expansions of the form 

d x A n = RA n + A n S 

so that we could prove that 

d x det(,4 n ) = Tr(i? + S) det(A n ). 

This search also proved fruitless, and in retrospect Lemma [3] suggests that it will be very difficult to 
find R and S because 

/ \ K 

(x -3)(2x -3){4x -3)Tr(R + S) = [8nx 2 - 6(5n + l)x - 3(9n 2 + 13n + 8) + 2(4n + 3) — - 

and experiments suggest that K n and H n are often relatively prime. 

So you have to feed some carefully selected facts into d x A ri to get a useful expansion of d x H n . 
Initially the first order differential equation 

x(x — l)d x a n (x) — (n(x — 3) — 2)a n {x) — (3n + 2)a ra (0) = 

looks like the fact needed. However this does not work because we never figured out what to do with 
the evaluation of the term 

Tr(A~ l [(Zi + 3j + 2)a i+j (0)] < i>:j < n ). 

So at this point we started searching for matrices that behaved well after being evaluated by the 
operator 

Tr^ 1 *) (165) 

By "behaves well" we mean that the calculations return a single determinant, or at most a linear 
combination of a just a few determinants, thus avoiding the expansion of H n as a sum of a large 
number of determinants. We quickly learned about the first two columns in Table [1] in Section [5j But 
this was not sufficient because d x A n cannot be expressed as a sum of the matrices defined in these 
two columns. 

At this point we recalled some results in the literature [HE] that related derivatives to convolutions. 
We realized that what this work told us was that the operator in (I165p works out very nicely on [cj+j+i], 
[cj +J ], [ci + j_i], i.e. the trace came out to be a single determinant as already shown in Tabled) 

This gave us enough tools to find our first identity, given as (jlip in Lemma[TJ Using this identity we 
have enough to express the first and second derivatives of H n as linear combinations of small numbers 
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of determinants. In particular in the examples of this paper, the first derivative of H n is expressed as 
a linear combination of H n and K n . The second derivative is expressed as a linear combination of H n , 
K n , M n and N n . This process can easily be continued to higher derivatives though the computations 
get messier at each derivative. 

If the first identity (Lemma [1]) is weakened to include terms like a n +2 and c n +2 the above process 
still works. The only difference is that the derivatives have more distinct determinant summands. The 
point of the above discussion is that we cannot handle any (/3, ct)-case in which there does not exist 
a first identity of this kind. So far there are a dozen cases, some not of the (/?, a) type, which have a 
first identity that we can handle, and we mention (3,0), (3,2), and (2,2) as examples of these. 

The role of the second (Lemma|2| and third identities ( I102H is to prove linear relationships between 
the determinants that are generated by the above process of differentiating H n . 

Of these three identities, the proof is most sensitive to the form of the second identity. In the 
(3, 2)-case, the second identity involves a n +3 terms. While this version of the second identity can still 
be used to prove linear relationships between determinants, there are not enough linear relationships, 
and a proof of the fourth order differential equation for the (3, 2)-case does not seem possible with the 
tools of this paper. 

Another difficulty that needs to be overcome is the problem of finding and proving the third 
identity. Our original approach was to guess the form of the third identity and thus reduce the 
problem to binomial identity proving techniques. This works quite well for the (2, l)-case where the 
third identity has an explicit form: 

w ( (n+l+ i\ (n+l+i\ ( (n + l + i\ (n+l + i\\ \ 

X> 1 > , ( 2 ( n + 1 _ i ) + ( n + 2 _ i ) + ( 2 ( n . t J + („ + i- i jj-j»- + -w=° cm 

for < m < n. 

Though we were eventually able to guess the third identity in the (3, l)-case, its form was hardly 
enlightening. In addition, the identity was so complex that even the job of applying automated tools 
to prove it would be a major undertaking. So the whole process had arrived at an impasse. 

The way forward was a sidestep via an existence theorem (Section H]). The third identity (once 
proved) said that the binomial identities we wished to prove actually existed, and were unique, and 
you only had to know a few of their components explicitly to show that the right hand side of (|96p 
vanished. The third identity for the (3, l)-case is (|102p . This identity is proved via the generating 
function of the a n (x) and from the explicit linear relationships governing the determinants H n (x), 
K n (x), N n (x) at special values of x. 
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Of special mention is the form of the generating function given in Lemma [5) It is important to 
know that this is not the form that emerges in the straightforward derivation of the generating func- 
tion, as was done in Lemma EJ We were surprised that we could not use this form directly to prove 
the third identity. It took a long time to discover the form in Lemma O This latter form worked, and 
this concluded the proof of the (3, l)-case. 

To sum up, a lot of things have to fall right in place perfectly to evaluate a H^ ,a \x). 

8 Additional results on Hankel determinants 

The proof technique presented here is applicable to Hankel determinants of polynomials other than the 
a}@> a \ Along with the (3, 0)-case, we give the necessary ingredients, i.e. the three identities required, 
for the proof of the differential equation satisfied for a few of these but omit the proofs of the theorems 
and the construction of the explicit power series solutions. 

Example 5 We start with the (3, 0)-case 
First identity for the (3, 0)-case: 



3(x — 3)x(4x — 3)d x a 



•n 



(4(2x - 3)n + 2(2x - 5))a„+i 



+(27(2x - 3)n + 3(4x 2 - 3x - 9))a, 
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(x - l)c n+ i + 27{x - l)Cn = 



Second identity for the (3, 0)-case: 



(4(2x - 3) 2 (5x - 3)n + 2(2x - 3) (5a - 3)(6x - ll))o n+2 

-(81(8a; 3 - 24x 2 + 27a; - 9)n + 18(37x 3 - 123x 2 + 153a; - 54))a n+1 

+(729x 3 ra + 486x 3 )a„ + 4(x - l)(2x - 3)(5x - 3)c n+2 

-3(40x 4 - 30x 3 - 207a; 2 + 270x - 81)c n+ i + 162x 2 (5x 2 - 15x + 9)c n = 



The generating function of the for the third identity for the (3, 0)-case: 

-(2x-3)t-3x 




Theorem 6 The Hankel determinant for the (3, 0)-case satisfies the differential equation 



{x - 3)(2x - 3)(5x - 3)d 2 x y - 2(10(n - l)x 2 - 9(3n - 4)x - 9(n + h))d x y 



+n(10(n - l)x - 3(n - 7))y = . 
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In the (3,0)-case, in addition to the already known product form at x = given in Q, at x = 3 
and surprisingly also at x = | the determinant is given by a simple product. We omit the proofs but 
record these evaluations below: 

(3)0) (3n)!(3n + 2)! -^ 3(6i-5)!(2i)!(2i-l) 
Hn (3)= 2(n!2) 11 (4i + l)!(4i-l)! ' (16?) 

" 27(6i-5)!(3i-l)(3i-2)(2i-l)! , . 

H * • ' * = 2 l 2 ( 4i- l).(4i- 3 )(4i- 4). ' (168 > 



Example 6 Next, take 

m^)=e nzri^ (i69) 



^2>k — 2m} 



m=0 N 

The first identity for the polynomials in (|169p : 



x(2x - 9)(4x + 9)4a„ - (36n + 30)a n+ i + (243n + 8x 2 + I8x + 81)a n 
-12c+i - (8x 2 - 36x - 81)c + 27x(2x - 9)c„_i = 

The second identity for the polynomials in (11691) : 

(36(2x + 3)n + 66(2x + 3))a n+2 - ((32x 3 + 486x + 729)n + 12(4x 3 + 4x 2 + 54x + 81))a n+ i 
+(216x 3 n + 108x 2 (2x + 3))a n + 12(2x + 3)c n+2 + (16x 3 - 72x 2 - 378x - 243)c n+ i 
+2x(8x 3 - 90x 2 + 243x + 729)c n - 54x 3 (2x - 9)c n _i = 

The generating function of the a k for the third identity for the polynomials in (|169|) : 

,, v 3t + 2x 

^' V ' ~ (x(4x - 9)y - 6)t + 9 + 2x - 6x 2 y 

where t 3 y = t — 1. 

Theorem 7 TTie Hankel determinant of the polynomials in \169\) satisfies the differential equation 

(2x + 3)(2x - 9)d 2 x y + 4(2(n + 2)x - 9(n + l))d x y - Y2n(n + l)y = . 
Example 7 What we refer to as the "aex"-case (for exceptional) is the Hankel determinant where 

a k (x) = -py £ ™) (m + l)(m + 2)x m (170) 



m=0 

The first identity for the aex-case: 



x(x — 3)(7x — 3)d x a n — (4(x — l)n + 14(x — l))a n+ i 

+(27(x - l)n + 3(x - l)(llx + 3))a„ - 6(x - l) 2 c n + 6x 3 c„_i = 
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The second identity for the aex-case: 

(4(x - l) 2 (3x - l)n + 18(x - l) 2 (3x - l))a n+2 + ((-113x 3 + 189x 2 - 135x + 27)n 
-4(30x 4 - 19x 3 - 21x 2 + 39x - 9))a„+i + (216x 3 n + 12x 3 (7x 2 - 2x + 15))a n 
+6(x - l) 3 (3x - l)c n+ i - 6x 2 (x - l)(Wx 2 - I7x + 9)c n + 6x 5 (7x - 9)c n _i = 

The generating function of the for the third identity for the aex-case: 

2(x - l) 2 r - 2(x - l) 2 - 2(3x - 1) 



f(x,y) 



(— 2yx 2 (x — 3)r + 2yx 2 (x — 3) + 6x 2 y — 3x + 1 
where 

r = T m y k . 

Theorem 8 The Hankel determinant of the polynomials in satisfies the differential equation 

(3x - l)(x - l)(x - 3)d 2 x y - 2(3nx 2 - 8nx - 3(n + 4))d x y + 3n(n + l)(x - l)y = . 

In the the determinant is given by a simple product for X rj CIS shown below: 

3 J r 2(6i + 7)!(2i + l)! 
(7) ~ JLi 7(4i + 5)!(4i + 3)! 

Example 8 We can aslo evaluate the Hankel determinants as an almost product for 

m=0 \ / 

These polynomials are related to the polynomials ajjf' ^ (x) of Example [5] by the simple transfor- 
mation 

m=0 \ / m=0 \ / 

Therefore 

Theorem 9 The Hankel determinant of the polynomials in ( 17$fy satisfies the differential equation 



(x - 2) (2a; - l)(5x + 2)4*/ - 2(10(n - l)z 2 - (7n - 16)x - 26n - 19)^2/ + n(10(n - 1> + 7n + ll)y = . 
Example 9 Finally we have an alternate evaluation of Hn'^ix) at x = 1. 
Theorem 10 Suppose afc(x) is defined as in ^ and H n (x) = det[aj+j(x)]o<ij<n- Then 

fr (6z + 4)!(2» + l)! " (-l)^!(4n + 3)!!(3n + i + 2)!(x-ir 

Li 2(4i + 2)!(4i + 3)! ^ (3n + 2)!i!(n - i)!(4n + 2« + 3)!! ' 1 j 
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As in Corollary U taking x = and x = 1 in (|174p we obtain the following alternate evaluations. 
Corollary 9 

(6i + 4)!(2i + l)! A n\{An + 3)!!(3n + i + 2)! 



det 



det 



3(* + j) + i 
i + j 

3(t + i) + 2 
i + j 



0<2j'<n 



n 



■E 



2(4i + 2)!(4i + 3)! ^ (3n + 2)H!(n - j)!(4n + 2« + 3)!! 



n 



(6i + 4)!(2i + l)! 



o< W <» i=1 2(4 l + 2)!(4 l + 3)! 

The first evaluation in Corollary [9] is the expression given in ([6]) for il^ . It is special as there are 
no cancellations on the right. The second evaluation is identical in form to the known product in (|3|) 
for^ 2) . 



There are other variants on the polynomials of' which experiments suggest satisfy differential 
equations. A particularly unusual example are the Hankel determinants for 




These determinants a third order equation, but the coefficients are very large and not round, making 
it hard to guess what they are. 

Going back to the Hankel determinants Hn , some of these are governed by a second order 
differential equation, such as (3,0), (3,1), and (2,1). There are also non-(/3, a)-cases governed by a 
second order differential equation such as the aex-case in Example [71 

In forthcoming work, we plan to develop the methods of this paper to investigate the Hankel 
determinants related to other families of polynomials. In the family (2, r), for example, the three cases 
r = 0,1,2 have product evaluations, but for r > 3, the evaluations are no longer products but only 
almost products. We also encounter new phenomena in these cases, such as higher order differential 
equations and case-splitting into residue classes. 

We remark that experimentally we know that both (2,3) and (3,2)-cases satisfy fourth order 
differential equations. For the (2, 3)-case, the polynomial in front of the fourth derivative has degree 
11. 



Acknowledgments: We would like to thank Julius Borcea, Ira Gessel, Christian Krattenthaler, Boris 
Shapiro, and Doron Zeilberger. 
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9 Appendix I: The generating function of the 

In Lemma 13 we give a closed form of the generating function f(x,y) defined in (|104p . 
Lemma 9 Suppose afc(x) is as defined in tlO\) and f is as in |ift^| j. Then 

/(X ' y)= (/3 + (l-^)(l-^) (175) 

where 

t< 3 y = t-l. (176) 
Proof Changing the order of summation and rearranging 

f( x ,y) = Y,* m y m j:( Pn + {(3 ~ 1)m + a )y n 

m>0 n>0 \ U J 

It is known [9] that 

g( „ > = ?T(T3fflI < 177 > 

where i satisfies ()176j) . For our generating function, (3 is the same but a in (|177p is replaced by 
(/3 — l)m + a. Using these parameters we obtain 

t {/3-l)m+a+l t a+l 

f(X,y)= So Xmyn 'P+O^Wt = ((3 + (1 - P)t)(l - xyt^) (178) 
where t satisfies ()176p . • 

Note that using the Lagrange inversion formula, t can be expanded as 
Proof of Lemma [5} 

Proof The expression (|178p for / with (3 = 3 and a = 1 can be rewritten in the form (|105p since 

i 2 4x + 2t - 6 

(3 - 2t)(l - xyi 2 ) ~ (x - 3)(4x - 3)y(2t - 3) - (x - l)(27y - 4) 

is equal to 

(t - 1 - t 3 y)(2tx + 6x - 9) 



(3 - 2t)(l - xyt 2 )(4tx 2 y - 6x 2 y - Ihtxy + 9xy + 9ty + 2x - 2) 
and the numerator has t — 1 — t 3 y as factor, but this is zero by (|176p . This gives the form of the 
generating function / for the (3, l)-case that is claimed in (I105p . • 
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Proof of Lemma [8} 

Proof The proof is similar to the proof of Lemma [5] once we observe that the expression for the 
generating function 

/('■»)- p-tia-M) (180> 

from Lemma [9] can be written as (|143|) since 

t 2 t 2t(l - t + t 2 y) 



(2-t)(l-xyt) (x - 2)yt + 1 - 2xy (t - 2)(1 - xyt)(l - 2yt - 2xy - xyt) 
and the right hand side vanishes since in this case (3 = 2 and 

l-t + t 2 y = (181) 

by Lemma [H • 
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10 Appendix II: Pairs of identities for the (3,1) and (2, l)-cases. 

Proofs of the first two identities in the (3, l)-case are as follows: 
Proof of Lemma [TJ 

Proof We make use of the generating function / = f(x, y) of the a^s in the form given by (I178p in 
Appendix I. Passing to the generating functions, (jlip is equivalent to the functional identity 

f — 1 f — 1 

(x - 3)(2x - 3)(4x - 3)d x f - 4yd v - 6- + (8x 2 - 18x + 36)/ (182) 

y y 

+27yd y f - 4(2x 2 - 6x + 3)/ 2 + 27(2x 2 - 6x + 3)yf = 



dyt = . (183) 



Using identity (|176p . 

f 
1 - 3yV 

Substituting this expression in the computation of d y f, the left hand side of (|182p can be simplified as 

2 ^ V ~ t ^ (2at 2 + 12x 2 yi 2 - 90xyt 2 + 54?/t 2 + 2t 2 + 8x 2 t 
(1 - 3yi 2 ) V 

+xt - 36x 2 yt + 135xyt - 8lyt - I5t - 8x 2 - 3x + 9 
which vanishes since 

t 3 y - t + 1 = 0by • 
Proof of Lemma [2} 

Proof We again use the generating function / = f(x, y) of the a^'s in the form given by (|178p . The 

identity (fl~3|) is equivalent to 

a , u , / -l-(4 + x)j/ ^ /-l-(4 + a)y 
8(x - l)ycL 5 h 20(x - 1) = 

y 2 y 2 

f — 1 f — 1 
-2(-81 + 135x - 72x 2 + lex 3 )^- 2(-117 + 180x - 92x 2 + 24x 3 )- 



f 2 _ I 

+27(2x - 3fyd y f + 54(2x - 3)(3 - 4x + 2x 2 )/ + 8(x - 1)(3 - 6x + 2x 2 ))—— (184) 



„2 11/|„3 , o„,4\ j-2 o7~/o™ oVn io„ i o™2\„,^2 



+2(81 - 297x + 324x 2 - 114x d + 8x*)f z - 27x(2x - 3)(9 - 12x + 2x')yf z = 

Using the expression in (|183p for d y t in the calculation of d y f, the left hand side of (|184p can be 
simplified as 

p-w^xl -**')» ( - 72(21,214 + 16tV * + 56W - 72>x ' - 1084 w 

-Wtx 3 - A8t 2 yx 3 - 156tyx 3 + 252yx 3 + 16x 3 - 4t 2 x 2 + 1134t 2 y 2 x 2 + 14tx 2 
+12t 2 yx 2 + 252tyx 2 - 486yx 2 - 10x 2 - 1701t 2 y 2 x + 32te + 18i 2 yx 
-i05tyx + 567yx - 24x + 4t 2 + 729t 2 y 2 - 30t - 54t 2 y + 243ty - 243y + 18 

which vanishes since t 3 y-t + 1 = by (JESJ. 
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Proofs of the first two identities in the (2, l)-case are as follows: 
Proof of Lemma [6} 

Proof Passing to the generating functions in (|130p . we need to prove the identity 

2x(x - 2)d x f - d y f + Ayd y f + 2{x + 1)/ - (x - l)/ 2 + 4(x - l)yf 2 = 0. 



Using the expression 



e- 



,y = — (185) 

in the calculation of d y f, the left hand side of the identity we want to prove can be simplified as 

2t 2 (t 2 y - t + 1)(2 - t + 2x + tty - 6txy) 
(2-t) 2 (l-2ty)(l-tey) 2 

which vanishes by (|18ip . • 
The proof of the second identity for the (2, l)-case is as follows: 

Proof of Lemma O 

Proof Again passing to the generating functions in (|13ip . we need to prove the identity 

f-l-{x + 3)y f - 1 - (x + 3)y f - 1 2 f-1 
xyd v h 2(x + 1) 5 2x(x + 2)yd v 2(2x z + 3x + 4) 

yZ yZ y y 

f 2 _ 1 

+ 8x 2 yd v f + 12x 2 f + (x - l)(x - 2)- 4(x - l)(x - 2)f 2 = 

y 

Using the expression for / and the expression for d y t in (|185p . the left hand side of this expression can 
be simplified to 

2(t 2 y-t + l) 



(2-t) 2 y 2 (l-2ty)(l-txyy 



■ (t 2 - 4 + Sty - 4t 2 y - 2t 3 y + 8txy - 2t 2 xy + t 3 xy 



-2t 2 x 2 y + 8t 3 y 2 - 16t 2 xy 2 - At 3 xy 2 + 8t 2 x 2 y 2 + 4t 3 x 2 y 2 - 4t 3 x 3 y 3 



which again vanishes by (|18ip . 
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11 Appendix III: On the degree of a class of Hankel Determinants 



Theorem 11 Let po,pi, ■ ■ ■ ,p n and qo, q\, . . . , q n be integer sequences. Further let 7 be real and 
Q!q, ai, . . . , a n and (3q, (3i, . . . , (3 n be sequences of real numbers. Then the determinant 



det 



E 



0<m<pi+qj 



K Pi + qj-m , 



(186) 



0<j,,j<n 



as a 



polynomial in x has degree < maxjmaxpj + max^j — n,0}. 



Proof We note the convention that the empty sum is zero, and the binomial coefficient in (I186p is 
interpreted via the gamma function. The proof of the theorem is by induction. At each stage of the 
induction we make use of the following lemma. 



Lemma 10 Let Po,Pi, ■ ■ ■ ,Pn an d Q01 Qi, ■ ■ ■ > Qn be two integer sequences and let ao, a>i,a,2> ■ ■ ■ be an 
infinite sequence of real numbers. Then the determinant 



det 



E 

0<m<pi+qj 



a Pi+qj—m x 



0<j,j'<n 

as a polynomial in x has degree < maxjmaxpj + maxqj — n, 0}. 

Proof By rearranging the indices, we can assume that the sequences po,pi, ■ ■ ■ ,p n and qo, q%, . . . , q n 
are nondecreasing without changing the conclusion of the lemma. Then 



det 



E a Pi+qj-m x 
0<m<pi+qj 



y 



T,Pi-T,ii .det 



0<i,j'<n 

where xy = 1. Therefore the analysis can focus on the degree of 



E 

0<m<pi+qj 



0<i,j<n 



det 



E 

0<m<pi+qj 



0<i ,j<n 

as a polynomial in y. By elementary row operations, we see that this is equal to the following 
determinant 

E ^ if' « ! 



det 



0<m<po+qj 

E a mD m otherwise 

0<m; pi-j+qj+lKmKpi+qj 



0<i,j<n 



Whenever p n + q n > n, the highest power of y in this determinant is no more than 

max{Efi +E^'°J' 
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and the lowest power of y in this determinant is at least 

n— 1 n— 1 

max{n + ^ + ^ <?j, 0} 
i=o j=0 

Now if we multiply this polynomial by y-Yv^Yn obtaining 



y 



det 



Er, 
a m y 

0<m<pi+qj 



0<i,j<n 

we get a polynomial in y . The degree of this polynomial is no more than 



max{p n + q n - n,0}. 

Converting back to x and remembering that pi and qj are nondecreasing gives the desired result. 
Now we are ready to start the proof of Theorem [TTJ Let 



CX ; 



With this change of variable, 



«i + (3j + 7m' 
.Pi + qj-m , 



r a[ + P'j - j(pi + qj - mf 
Pi + qj-m 



The following identity holds: 

\ Pi + qj -m , 



E E 

r>0 s=0 



K Pi — r + qj — s — m. 



This allows us to express the determinant as 



det 



EE HP) E ( 

r>0s>0 \ / \ / 0<m<pj— r+gj— s V ' 



7(Pi + ??-"») 



0<i,j<n 



This determinant has an expansion as a sum of terms of the form 



11 ; 11 



'0: 



.det 



E 

0<m<pi— n+qj— Sj 



-jPi + -jqj + jm \ m 
a>i -n + qj - Sj-mJ 



0<i,j<n 



over collections of nonnegative integers tq, r^, . . . , v n and So? s\, . . . , s n . 

If all the r, and Sj in any such collection are zero, then the corresponding determinant is as 
described in Lemma [lOl and therefore has the appropriate degree as a polynomial in x. If any of the 
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rj or Sj are nonzero, then the determinant has the same form as that of Theorem [TT] except that pi 
has been replaced with p' i = p% — r« and q'j = qj — Sj. This is the induction process. 
Note that the induction process terminates, since eventually 

max{maxj)' + max q'j — n, 0} = 

and in fact the resulting determinants in (|186p become zero. • 

Corollary 10 The + x (n+1) Hankel determinant Hn'°^ defined by the sequence of polynomials 
in {Di is of degree at most n. 

Proof Hn ,a - 1 is of the type described in Theorem [TT1 where p% = q% = i, on = (3i + a, Pi = (3i, and 
7 = —1. • 
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